A Guided Walk Down Wall Street: an Introduction to Econophysics 
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This article contains the lecture notes for the short course "Introduction to Econophysics," de- 
hvered at the II Brazilian School on Statistical Mechanics, held in Sao Carlos, Brazil, in February 
2004. The main goal of the present notes is twofold: i) to provide a brief introduction to the problem 
of pricing financial derivatives in continuous time; and ii) to review some of the related problems to 
which physicists have made relevant contributions in recent years. 



I. INTRODUCTION 

This article comprises the set of notes for the short 
course "Introduction to Econophysics," dehvered at the 
II Brazihan School on Statistical Mechanics, held at the 
University of Sao Paulo, in Sao Carlos, SP, Brazil, in 
February 2004. The course consisted of five lectures and 
was aimed at physics graduate students with no previous 
exposition to the subject. 

The main goal of the course was twofold: i) to provide a 
brief introduction to the basic models for pricing financial 
derivatives; and ii) to review some of the related problems 
in Finance to which physicists have made significant con- 
tributions over the last decade. The recent body of work 
done by physicists and others have produced convincing 
evidences that the standard model of Finance (see be- 
low) is not fully capable of describing real markets, and 
hence new ideas and models are called for, some of which 
have come straight from Physics. In selecting some of 
the more recent work done by physicists to discuss here, 
I have tried to restrict myself to problems that may have 
a direct bear on models for pricing derivatives. And even 
in such cases only a brief overview of the problems is 
given. It should then be emphasized that these notes are 
not intended as a review article on Econophysics, which 
is nowadays a broad interdisciplinary area, but rather 
as a pedagogical introduction to the mathematics (and 
physics?) of financial derivatives. Hence no attempt has 
been made to provide a comprehensive list of references. 

No claim of originality is made here regarding the 
contents of the present notes. Indeed, the basic the- 
ory of financial derivatives can now be found in numer- 
ous textbooks, written at a different mathematical lev- 
els and aimingat specific (or mixed) audiences, such as 
economists [ p, 13. IE 1 11 , apphed mathematicians 0,13013; 
physicists 0, Hfl HHi etc. (Here I have listed only the 
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texts that were most often consulted while writing these 
notes.) Nevertheless, some aspects of presentation given 
here have not, to my knowledge, appeared before. An 
example is the analogy between market efficiency and a 
certain symmetry principle that is put forward in Sec. V. 
Similarly, the discussion of some of the more recent re- 
search problems is based on the already published litera- 
ture. An exception is Fig. 1 121 which contains unpublished 
results obtained by R. L. Costa and myself. 

The present notes are organized as follows. Section II 
gives some basic notions of Finance, intended to intro- 
duce the terminology as well as the main problems that 
I shall be considering. In Sec. Ill, I discuss the Brow- 
nian motion, under a more formal viewpoint than most 
Physics graduate students are perhaps familiar with, and 
then develop the so-called Ito stochastic calculus. Section 
IV contains what is the raison d'etre of the present notes, 
the Black-Scholes model for pricing financial derivatives. 
In Sec. V, the martingale approach for pricing deriva- 
tives is introduced. In particular, I recast the notions of 
market efficiency and no-arbitrage as a 'symmetry prin- 
ciple' and its associated 'conservation law.' Sections VI 
and VII discuss two possible ways in which real markets 
may deviate from the standard Black-Scholes model. The 
first of such possibilities is that financial asset prices have 
non-Gaussian distributions (Sec. VI), while the second 
one concerns the presence of long-range correlations or 
memory effects in financial data (Sec. VII). Conclusions 
are presented in Sec. VIII. For completeness, I give in 
Appendix A the formal definitions of probability space, 
random variables, and stochastic processes. 



II. BASIC NOTIONS OF FINANCE 

A. Riskless and risky financial assets 

Suppose you deposit at time t — an amount of R$ 1 
into a bank account that pays an interest rate r. Then 
over time the amount of money you have in the bank, let 



2 



us call it B{t), will increase at a rate 



Solving this equation subject to the initial condition 
B(0) = 1 yields 

B{t)^e^K (2) 

A bank account is an example of a riskless financial 
assets, since you are guaranteed to receive a known (usu- 
ally fixed) interest rate r, regardless of the market situ- 
ation. Roughly speaking, the way banks operate is that 
they borrow from people who have money to 'spare', but 
are not willing to take risks, and lend (at higher interest 
rates) to people who 'need' money, say, to invest in some 
risky enterprise. By diversifying their lending, banks can 
reduce their overall risk, so that even if some of these 
loans turn bad they can still meet their obligations to 
the investors from whom they borrowed. 

Governments and private companies can also borrow 
money from investors by issuing bonds. Like a bank ac- 
count, a bond pays a (fixed or floating) interest rate on 
a regular basis, the main difference being that the repay- 
ment of the loan occurs only at a specified time, called the 
bond maturity. Another difference is that bonds are not 
strictly risk-free assets because there is always a chance 
that the bond issuer may default on interest payments 
or (worse) on the principal. However, since governments 
have a much lower risk to default than corporations, cer- 
tain government bonds can be considered to be risk free. 

A company can also raise capital by issuing stocks or 
shares. Basically, a stock represents the ownership of 
a small piece of the company. By selling many such 
'small pieces', a company can raise capital at lower costs 
than if it were to borrow from a bank. As will be dis- 
cussed shortly, stocks are risky financial assets because 
their prices are subjected to unpredictable fluctuations. 
In fact, this is what makes stocks attractive to aggressive 
investors who seek to profit from the price fluctuations 
by pursuing the old advice to "buy low and sell high." 

The buying and selling of stocks are usually done 
in organized exchanges, such as, the New York Stock 
Exchange (NYSE) and the Sao Paulo Stock Exchange 
(BOVESPA). Most stock exchanges have indexes that 
represent some sort of average behavior of the corre- 
sponding market. Each index has its own methodology. 
For example, the Dow Jones Industrial Average of the 
NYSE, which is arguably the most famous stock index, 
corresponds to an average over 30 industrial companies. 
The Ibovespa index of the Sao Paulo Stock Exchange, 
in contrast, represents the present value of a hypotheti- 
cal portfolio made up of the stocks that altogether corre- 
spond to 80% of the trading volume. Another well known 
stock index is the Standard & Poor's 500 (S&P500) In- 
dex calculated on the basis of data about 500 companies 
listed on the NYSE. [Many other risky financial assets, 
such as, currency exchange rates, interest rates, and com- 
modities (precious metals, oil, grains, etc), are traded on 
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FIG. 1: Daily closing values of the deflated Ibovespa index in 
the period 1968-2003. 



organized markets but these will not be discussed any 
further in the present notes.] 

B. The random nature of stock prices 

Since a stock represents a 'small piece' of a company, 
the stock price should somehow reflect the overall value 
(net worth) of this company. However, the present value 
of a firm depends not only on the firm's current situ- 
ation but also on its future performance. So here one 
sees already the basic problem in pricing risky financial 
assets: we are trying to predict the future on the ba- 
sis of present information. Thus, if a new information is 
revealed that might in one way or another affect the com- 
pany's future performance, then the stock price will vary 
accordingly. It should therefore be clear from this simple 
discussion that the future price of a stock will always be 
subjected to a certain degree of uncertainty. This is re- 
flected in the typical 'erratic behavior' that stock prices 
show when graphed as a function of time. An example 
of such a graph is shown in Fig. ^ for the case of the 
Ibovespa stock index. 

Although stock prices may vary in a rather unpre- 
dictable way, this does not mean that they cannot be 
modeled. It says only that they should be described in a 
probabilistic fashion. To make the argument a little more 
precise, let S be the price of a given stock and suppose 
we want to write an equation analogous to for the 
rate of increase of S: 

^ = ms, (3) 

where R{t) represents the 'rate of return' of the stock. 
The question then is: what is R{t)7 From our previous 
discussion, it is reasonable to expect that R{t) could be 
separated into two components: i) a predictable mean 
rate of return, to be denoted by /i, and ii) a fluctuat- 
ing ('noisy') term ^(t), responsible for the randomness 
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or uncertainty in the stock price. Thus, after writing 
R{t) = ^ + ^(t) in © we have 

^ = [^, + ^{t)]s. (4) 

Now, one of the best models for 'noise' is, of course, the 
white noise, so it should not come as a surprise to a 
physicist that Brownian motion and white noise play an 
important role in finance, as will be discussed in detail 
shortly. 

C. Options and derivatives 

Besides the primary financial assets already mentioned 
(stocks, commodities, exchange rate, etc), many other 
financial instruments, such as options and futures con- 
tracts, are traded on organized markets (exchanges). 
These securities are generically called derivatives, be- 
cause they derive their value from the price of some pri- 
mary underlying asset. Derivatives are also sometimes 
referred to as contingent claims, since their values are 
contingent on the evolution of the underlying asset. In 
the present notes, I will discuss only one of the most basic 
derivatives, namely, options. 

An option is a contract that gives its holder the right, 
but not the obligation, to buy or sell a certain asset for a 
specified price at some future time. The other part of the 
contract, the option underwriter, is obliged to sell or buy 
the asset at the specified price. The right to buy (sell) is 
called a call (put) option. If the option can only be ex- 
ercised at the future date specified in the contract, then 
it is said to be a European option. American options, on 
the other hand, can be exercised at any time up to matu- 
rity. (For pedagogical reasons, only European derivatives 
will be considered here.) To establish some notation let 
us give a formal definition of a European option. 

Definition 1 A European call option with exercise price 
( or strike price ) K and maturity ( or expiration date ) T 
on the underlying asset S is a contract that gives the 
holder the right to buy the underlying asset for the price 
K at time T . 

A European put option is the same as above, the only 
difference being that it gives the holder the right to sell 
the underlying asset for the exercise price at the expira- 
tion date. 

If at the expiration date T the stock price St is above 
the strike price K , the holder of a call option will exercise 
his right to buy the stock from the underwriter at price K 
and sell it in the market at the spot price St, pocketing 
the difference St — K. On the other hand, if at expiration 
the price St closes below K then the call option becomes 
worthless (since it would be cheaper to buy the stock in 
the market). The payoff oi a call option at maturity is 
therefore given by 

payoff call = niax(S'T -K,{)). (5) 
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FIG. 2: Value of a call option at the expiration date (thick 
line) and before expiration (thin line). 
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FIG. 3: Value of a put option at the expiration date (thick 
line) and before expiration (thin line). 



The payoff diagram of a call option is illustrated by the 
thick line in Fig. [21 In this figure the thin line represents 
the price of the call option at an arbitrary time t < T 
before expiration. (The otpion price before expiration is 
always greater than the payoff at expiration on account 
of the higher risks: the further way the expiration date, 
the greater the uncertainty regarding the stock price at 
expiration.) Similarly, the payoff function for a put op- 
tion is 

Payoffput = niax(X - 5t, 0), (6) 

which is shown as the thick line in Fig. |3| 

Because an option entitles the holder to a certain right 
it commands a premium. Said differently, since the un- 
derwriter has an obligation (while the holder has only 
rights) he will demand a payment, to be denoted by Co, 
from the holder in order to enter into such a contract. 
Thus, in the case of a call option, if the option is exer- 
cised the holder (underwriter) will make a profit (loss) 
given by max(5' — K,0) — Cq; otherwise, the holder (un- 
derwriter) will have lost (won) the amount Cq paid (re- 
ceived) for the option. And similarly for a put option. 
Note then that the holder and the underwriter of an op- 
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tion have opposite views regarding the direction of the 
market. For instance, the holder of a call option is bet- 
ting that the stock price will increase (past the exercise 
price), whereas the underwriter hopes for the opposite. 

Now, given that the holder and the underwriter have 
opposite views as to the direction of the market, how can 
they possibly agree on the price for the option? For if the 
holder (underwriter) suspects that the option is overval- 
ued (undervalued) he will walk away from the contract. 
The central problem in option pricing is therefore to de- 
termine the rational price price Cq that ensures that nei- 
ther part 'stands a better chance to win.' 

A solution to this problem (under certain assumptions) 
was given in 1973 in the now-famous papers by Black 
and Scholes and Merton , which won Scholes and 
Merton the Nobel prize in Economics in 1997. (Black had 
died meanwhile.) The history of options is however much 
longer. In fact, the first scientific study of options dates 
back to the work by the French mathematician Bachelier 
in 1900 who solved the option pricing problem above 
but under slightly wrong assumptions; see, e.g., (ll| for 
a detailed discussion of Bachelier's work. 

After an option (traded on exchange) is first underwrit- 
ten, it can subsequently be traded and hence its 'market 
price' will be determined by the usual bid-ask auction. It 
is nonetheless important to realize that investors in such 
highly specialized market need some basic pricing the- 
ory to rely on, otherwise investing in options would be 
a rather wild (and dangerous) game. Indeed, only after 
the appearance of the Black-Scholes model [and the es- 
tablishment of the first option exchange in Chicago also 
in 1973] have option markets thrived. One of the main 
objectives of the present notes is to explain the theo- 
retical framework, namely, the Black-Scholes model and 
some of its extensions, in which options and other deriva- 
tives are priced. I will therefore not say much about the 
practical aspects of trading with options. 



D. Hedging, speculation, and arbitrage 

Investors in derivative markets can be classified into 
three main categories: hedgers, speculators, and arbi- 
trageurs. 

Hedgers are interested in using derivatives to reduce 
the risk they already face in their portfolio. For example, 
suppose you own a stock and are afraid that its price 
might go down within the next months. One possible 
way to limit your risk is to sell the stock now and put 
the money in a bank account. But then you won't profit 
if the market goes up. A better hedging strategy would 
clearly be to buy a put option on the stock, so that you 
only have to sell the stock if it goes below a certain price, 
while getting to keep it if the price goes up. In this case 
an option works pretty much as an insurance: you pay 
a small price (the option premium Cq) to insure your 
holdings against possibly high losses. 



Speculators, in contrast to hedgers, seek to make profit 
by taking risks. They 'take a position' in the market, 
by betting that the price on a given financial asset will 
go either up or down. For instance, if you think that 
a certain stock will go up in the near future, you could 
"buy and hold" the stock in the hope of selling it later 
at a profit. But then there is the risk that the price goes 
down. A better strategy would thus be to buy a call 
option on the stock. This not only is far cheaper than 
buying the stock itself but also can yield a much higher 
return on your initial investment. (Why?) However, if 
the market does not move in the way you expected and 
the option expire worthless, you end up with a 100% 
loss. (That's why speculating with option is a very risky 
business.) 

Arbitrageurs seek to make a riskless profit by entering 
simultaneously into transactions in two or more markets, 
usually without having to make any initial commitment 
of money. The possibility of making a riskless profit, 
starting with no money at all, is called an arbitrage op- 
portunity or, simply, an arbitrage. A more formal defini- 
tion of arbitrage will be given later. For the time being, 
it suffices to give an example of how an arbitrage oppor- 
tunity may arise. 

But before going into this example, it is necessary first 
to discuss the notion of a short sell. 'Shorting' means 
selling an asset that one does not own. For example, 
if you place an order to your broker to short a stock, 
the broker will "borrow" a stock from somebody else's 
account, sell it in the market, and credit the proceeds into 
your account. When you then decide to close your short 
position (there usually is a limit on how long an asset 
can be held short), your broker will buy the stock in the 
market (taking the money from your account) and return 
it to its original owner. If in the meantime the stock 
prices decreased, the short sell brings a profit, otherwise 
the short seller incurs in a loss. This is why a short sell 
is usually done simultaneously with another operation 
to compensate for this risk (as in the arbitrage example 
below). It should also be noted, in passing, that buying 
the actual asset corresponds to taking a 'long position' 
on this asset. 

Let us now consider our hypothetical arbitrage exam- 
ple. Many Brazilian companies listed in the Sao Paulo 
Stock Exchange also have their stocks traded on the New 
York Stock Exchange in the form of the so-called Amer- 
ican Depository Receipt (ADR). Suppose then that a 
stock is quoted in Sao Paulo at R$ 100, with its ADR 
counterpart trading in New York at US$ 34, while the 
currency rate exchange is 1 USD = 2.90 BRL. Start- 
ing with no initial commitment, an arbitrageur could sell 
short N stocks in Sao Paulo and use the proceeds to buy 
N ADR's in New York (and later have them transferred 
to Sao Paulo to close his short position). The riskless 
profit in such operation would be R$ (100-2.90 x 34)iV = 
R$ 1.40 N . (In practice, the transaction costs would elim- 
inate the profit for all but large institutional investors 

Q.) 
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p= 1/2 



Sn = 57 




p= 1/2 



Si=65 



Si=53 



FIG. 4: One-step binomial model for a stock. 



Vo = Co-57A 




p= 1/2 



Vi = 8-65A 



Vi = -53A 



FIG. 6: Delta-hedging portfolio in the one-step binomial 
model. 



p= 1/2 




Ci=65-57 = 8 



p= 1/2 



Ci = 



FIG. 5: Option value in the one-step binomial model. 



Note, however, that such 'mispricing' cannot last long: 
buy orders in New York will force the ADR price up, 
while sell orders in Sao Paulo will have the opposite ef- 
fect on the stock price, so that an equilibrium price for 
both the ADR and the stock is soon reached, whereupon 
arbitrage will no longer be possible. In this sense, the 
actions of an arbitrageur are self-defeating, for they tend 
to destroy the very arbitrage opportunity he is acting 
upon — but before this happens a lot of money can be 
made. Since there are many people looking for such risk- 
less chances to make money, a well-functioning market 
should be free of arbitrage. This is the main idea behind 
the principle that in an efficient market there is no arbi- 
trage, which is commonly known as the "no-free-lunch" 
condition. 



E. The no-arbitrage principle in a (binomial) 
nutshell 



Here we shall consider a one-step binomial model to 
illustrate the principle of no-arbitrage and how it can be 
used to price derivatives. Suppose that today's price of 
an ordinary Petrobras stocks (PETR3 in their Bovespa 
acronym) is 5*0 — 57 BRL. Imagine then that in the next 
time-period, say, one month, the stock can either go up 
to Si = 65 with probability p or go down to Sf = 53 with 
probability q. For simplicity let us takep = q = 1/2. Our 



binomial model for the stock price dynamics is illustrated 
in Fig. 01 Note that in this case the stock mean rate of 
return, /i, is given by the expression: (1 -I- /Lt)S'o = £^[5"!], 
where £^[5*] denotes the expected value of S (see Sec. lIII Al 
for more on this notation). Using the values shown in 
Fig. ^ one then gets — 0.035 or fi — 3.5%. Let us 
also assume that the risk-free interest rate is r = 0.6% 
monthly. 

Consider next a call option on PETR3 with exercise 
price K = 57 and expiration in the next time period, 
i.e., T — 1. Referring to (|3J) and Fig. 21 one immediately 
sees that at expiration the possible values (payoffs) for 
this option in our binomial model are as shown in Fig. [31 
C" = 8 or Cf = with equal probability. The ques- 
tion then is to determine the 'rational' price Co that one 
should pay for the option. Below we will solve this prob- 
lem using two different but related methods. The idea 
here is to illustrate the main principles involved in op- 
tion pricing, which will be generalized later for the case 
of continuous time. 

First, we describe the so-called delta- hedging argu- 
ment. Consider a portfolio made up of one option C 
and a short position on A stocks, where A is to be deter- 
mined later, and let Vt denote the money value of such a 
portfolio at time t. We thus have 

Vt = Ct- ASu 

where the minus sign denotes that we have short sold 
A stocks (i.e., we 'owe' A stocks in the market). From 
Figs. 21 and |S1 one clearly sees that the possibles values 
for this portfolio in our one-step model are as illustrated 
in Fig. |B1 Let us now chose A such that the value Vi 
of the portfolio is the same in both 'market situations.' 
Referring to Fig.jSlone immediately finds 



Vi = Vi =^ 8 - A ■ 65 = -A ■ 5S 



A = l 
3 



Thus, by choosing A = 2/3 we have completely elimi- 
nated the risk from our portfolio, since in both (up or 
down) scenarios the portfolio has the same value Vi . But 
since this portfolio is riskless, its rate of return must 
be equal to the risk-free interest rate r, otherwise there 
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would be an arbitrage opportunity, as the following ar- 
gument shows. 

Let r' denote the portfolio rate of return, i.e., r' is the 
solution to the following equation 

{l + r')Vo^Vi. (7) 

If r' < r, then an arbitrageur should take a long position 
on (i.e., buy) the option and a short position on A stocks. 
To see why this is an arbitrage, let us go through the 
argument in detail. At time t = the arbitrageur's net 
cashflow would he Bq — |Vb| = A • — Co, which he 
should put in the bank so that in the next period he 
would have Bi = (1 + r)|Vo|. At time t = 1, he should 
then close his short position on A stocks, either exercising 
his option (up scenario) or buying A stocks directly in 
the market (down scenario). In either case, he would 
have to pay the same amount \Vi\ = (1 + r')|Vb| < -Bi, 
and hence would be left with a profit of Bi — \Vi\. On 
the other hand, if r' > r the arbitrageur should adopt 
the opposite strategy: go short on (i.e., underwrite) the 
option and long on A stocks (borrowing money from the 
bank to do so). 

We have thus shown that to avoid arbitrage we must 
have r' — r. This is indeed a very general principle that 
deserves to be stated in full: in a market free of arbitrage 
any riskless portfolio must yield the risk-free interest rate 
r. This no-arbitrage principle is at the core of the modern 
theory of pricing derivatives, and, as such, it will be used 
several times in these notes. 

Let us now return to our option pricing problem. Set- 
ting r' = r in Q and substituting the values of Vq and 
Vi given in Fig. El we obtain 

{l + r)[Co-ASo] = -ASf. (8) 

Inserting the values of r = 0.006, Sq = 57, Sf = 53, and 
A = 2/3 into the equation above, it then follows that the 
option price that rules out arbitrage is 

Co = 2.88. (9) 

It is instructive to derive the option price through a 
second method, namely, the martingale approach or risk- 
neutral valuation. To this end, we first note that from 
Fig. |S1 we see that the expected value of the option at 
expiration is i?[Ci] = ^8-1-^0 = 4. One could then 
think, not totally unreasonably, that the correct option 
price should be the expected payoff discounted to the 
present time with the risk-free interest rate. In this case 
one would get 

C-M^^^ = 3.98, 
" l + r 1.006 

which is quite different from the price found in Q. The 
faulty point of the argument above is that, while we used 
the risk-free rate r to discount the expected payoff £'[Ci], 
we have implicitly used the stock mean rate of return fi 
when calculating i?[Ci]. Using these two different rates 



leads to a wrong price, which would in turn give rise to 
an arbitrage opportunity. 

A way to avoid this arbitrage is to find fictitious prob- 
abilities Qu and Qd, with qu + qa = 1, such that the stock 
expected return calculated with these new probabilities 
would equal the risk- free rate r. That is, we must de- 
mand that 

So{l + r)=E'^ [Si] = g„ • S'l + qa ■ , (10) 

where iS'^Jx] denotes expected value with respect to the 
new probabilities qu and qd- Using the values for Si and 
Sf given in Fig. ^ we easily find that 

qu = 0.3618, qd = 0.6382. 

Under these probabilities, the expected value iJ'^JC'j^] of 
the option at maturity becomes E'^[Ci\ = 0.3618 x 8 -I- 
0.6382 X = 2.894, which discounted to the present time 
yields 

_ EQ[Ci] 2.894 

^° - TT^ - 1006 - 

thus recovering the same price found with the delta- 
hedging argument. 

Note that under the fictitious probability <7„ and qd, 
all financial assets (bank account, stock, and option) in 
our binomial model yield exactly the same riskless rate 
r. Probabilities that have this property of 'transforming' 
risky assets into seemingly risk-free ones are called an 
equivalent martingale measure. Martingale measures is a 
topic of great relevance in Finance, as will be discussed 
in more detail in Sec. IV. 

In closing this subsection, it should be noted that the 
one-step binomial model considered above can be easily 
generalized to a binomial tree with, say, N time steps. 
But for want of space this will not be done here. (I 
anticipare here, however, that Black-Scholes model to be 
considered later corresponds precisely to the continuous- 
time limit of the binomial multistep model.) It is perhaps 
also worth mentioning that binomial models are often 
used in practice to price exotic derivatives, for which no 
closed formula exists, since such models are rather easy 
to implement on the computer; see, e.g., for more 
details on binomial models. 



F. Put-Call parity 

In the previous subsection I only considered the price of 
a (European) call option, and the attentive reader might 
have wondered how can one determine the price of the 
corresponding put option. It turns out that there is a 
simple relationship between European put and call op- 
tions, so that from the price of one of them we can ob- 
tain the price of the other. To see this, form the following 
portfolio: i) buy one stock 5* and one put option P on 
this stock with strike price K and maturity T, and ii) 
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short one call option C with the same strike and matu- 
rity as the put option. The value of such portfolio would 
thus be 

V = S + P-C. (11) 

Now from and (|SJ), one immediately sees that at ex- 
piration we have P — C = K — S, so that the value of the 
above portfolio at time T becomes simply 

Vt = K. (12) 

Since this portfolio has a known (i.e., riskless) value at 
time t — T, it then follows from the no-arbitrage condi- 
tion that its value at any time < t < T must be given 
by 

V^Ke''-^^-*\ (13) 

where r is the risk-free interest rate. Inserting (|13|) into 
((mi immediately yields the so-called put-call parity re- 
lation: 

P = C - S + Ke-'''^^-*\ (14) 

III. BROWNIAN MOTION AND STOCHASTIC 
CALCULUS 

A. One-dimensional random walk 

Every physics graduate student is familiar, in one way 
or another, with the concept of a Brownian motion. The 
customary introduction [l5l | to this subject is through the 
notion of a random walk, in which the anecdotal drunk 
walks along a line taking at every time interval At one 
step of size I, either to the right or to the left with equal 
probability. The position, X{t), of the walker after a time 
t = NAt, where N is the number of steps taken, repre- 
sents a stochastic process. (See Appendix A for a formal 
definition of random variables and stochastic processes.) 
As is well known, the probability P{X{t) = x) for the 
walker to be found at a given position x = nl, where n 
is an integer , at given time t, is described by a binomial 
distribution 

Simply stated, the Brownian motion is the stochastic 
process that results by taking the random walk to the 
continuous limit: At —fO^l—^O^N—^ oo, n — > oo 
such that t — NAt and x = nl remain finite. (A more 
formal definition is given below.) Here, however, some 
caution with the limits must be taken to ensure that a 
finite probability density p(x, t) is obtained: one must 
take At ^ Q and Z ^ 0, such that P = a At, where a is 
a constant. In this case one obtains that p{x,t) is given 
by a Gaussian distribution p^ : 

V2T:a^t L 2cr^t J 

At this point let us establish some notation. Let X 
be a random variable with probability density function 



(pdf) given by p{x). [Following standard practice, we 
shall denote a random variable by capital letters, while 
the values it takes will be written in small letters]. The 
operator for expectation value will be denoted either as 
E[-] or < • >, that is, 

/oo 
f{x)p{x)dx, (16) 
-oo 

where f{x) is an arbitrary function. Although the 
angular-bracket notation for expectation value is pre- 
ferred by physicists, we shall often use the E notation 
which is more convenient for our purposes. 

A Gaussian or normal distribution with mean m and 
standard deviation a will be denoted by J\f{m, cr), whose 
pdf is 

Let us also recall that the (nonzero) moments of the 
Gaussian distribution are as follows 

E[X] = m, E[X^] = CT^ (18) 
EiX^""] = 1-3-5-... •(2n-l)cr2". (19) 

B. Brownian motion and white noise 

We have seen above that a ID Brownian motion can 
be thought of as the limit of a random walk after in- 
finitely many infinitesimal steps. This formulation was 
first given in 1900 by Bachelier who also showed the 
connection between Brownian motion and the diffusion 
equation (five years before Einstein's famous work on the 
subject It is thus telling that the first theory of 

Brownian motion was developed to model financial asset 
prices! A rigorous mathematical theory for the Brown- 
ian motion was constructed by Wiener [l7l| in 1923, after 
which the Brownian motion became also known as the 
Wiener process. 

Definition 2 The standard Brownian motion or Wiener 
process {W{t),t > 0} is a stochastic process with the fol- 
lowing properties: 

1. W{Q) = 0. 

2. The increments W{t) — W{s) are stationary and 
independent. 

3. For t > s, W{t) — W{s) has a normal distribution 

4-. The trajectories are continuous (i.e., "no jumps"). 

The stationarity condition implies that the pdf of 
W{t) — W{s), for t > s, depends only on the time differ- 
ence t — s. (For a more precise definition of stationary 
processes see Appendix A.) Now, it is not hard to con- 
vince oneself that conditions 2 and 3 imply that W{t) is 
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distributed according to A^(0, y/t) for t > 0. In partic- 
ular, we have = for all i > 0. Furthermore, 
one can easily show that the covariance of the Brownian 
motion is given by 

E[W{t)W{3)] = s, for t > s. 

It is also clear from the definition above that the Brow- 
nian motion is a Gaussian process (see Appendix A for 
the formal definition of Gaussian processes). Then, since 
a Gaussian process is fully characterized by its mean and 
covariance, we can give the following alternative defini- 
tion of the Brownian motion. 

Definition 3 The standard Brownian motion or Wiener 
process {W{t), t > 0} is a Gaussian process with 
E[W{t)] = and E[W{t)W{s)] = min(s,i). 

The Brownian motion has the important property of 
having bounded quadratic variation. To see what this 
means, consider a partition {ti}"^Q of the interval [0,t], 
where = io < < • ■ • < = For simplicity, let us 

take equally spaced time intervals: — ti_i = At = — . 

n 

The quadratic variation of W{t) on [0, t] is defined as 

n 

g„ = ^AWi2, (20) 



i=0 



where AW^ = W{t,)-W{t^.i). Since AWi is distributed 
according to 7V(0, y/Ki) we have that ElAW"^] = At, 
which implies that 



E[Qn] - t. 



(21) 



Furthermore, using the fact that the increments AWi are 
independent and recalling that the variance of the sum 
of independent variables is the sum of the variances, we 
get for the variance of Q„: 

n n 

var[Q„] ^ 5]var[AM/f] = {^[^^'] ^ (mW^)'} 



i=0 



j=0 



Y: [3(At)2 - [Atf] 



2t^ 



i=0 



where in the third equality we used (|19|l and the fact that 
AWi has distribution J\f{0, VAt). We thus see that 



var[Qn] — > 0, as n 
On the other hand, we have that 

^2 



var[Q„]=E (Qn-E[Qn])^ 



E 



(Qn-t)^ 



(22) 



(23) 



where in the last equality we have used ()21fl . Comparing 
(1^ and ^ then yields 



lim E 



{Qn - ty 



= 0. 



We have thus proven that Qn converges to t in the 
mean square sense. This fact suggests that AW'^ can be 
thought of as being of the order of At, meaning that as 
At the quantity AW'^ "resembles more and more" 
the deterministic quantity At. In terms of differentials, 
we write 

[dW]"^ = dt. (24) 
Alternatively, we could say that dW is of order \/dt: 



dW ^ 0{Vdt) 



(25) 



(I remark parenthetically that the boundedness of the 
quadratic variation of the Brownian motion should be 
contrasted with the fact that its total variation. An = 
J2i=o is unbounded, that is. An ^ cxd as n ^ oo, 

with probability 1; see 0.) 

Another important property of the Brownian motion 
W(t) is the fact that it is self-similar (or more exactly 
self-ajjine) in the following sense: 



W{at) = a^/^W{t), 



(26) 



for all a > 0. Here — means equality in the sense of prob- 
ability distribution, that is, the two processes W{at) and 
a^^'^W{t) have exactly the same finite-dimensional distri- 
butions p{xi, ti] ...,Xn,tn) for any choice of t^, z = 1, rt, 
and n > 1. Self-similarity means that any finite portion 
of a Brownian motion path when properly rescalcd is 
(statistically) indistinguishable from the whole path. For 
example, if we 'zoom in' in any given region (no matter 
how small) of a Brownian motion path, by rescaling the 
time axis by a factor of a and the vertical axis by a factor 
of -y/a, we obtain a curve similar (statistically speaking) 
to the original path. An example of this is shown in 
Fig- [3 In the language of fractals, we say that a trajec- 
tory of a Brownian motion is a fractal curve with fractal 
dimension D ^ 2. 

The self-similarity property implies that sample paths 
of a Brownian motion are nowhere differentiable (tech- 
nically, with probability 1). A formal proof of this fact, 
although not difficult, is beyond the scope of the present 
notes, so that here we shall content ourselves with the fol- 
lowing heuristic argument. Suppose we try to compute 
the derivative of W{t) in the usual sense, that is. 



lim 



AW 



dW 

dt At'-^o At 



lim 

At^Q 



W{t + At) - W{t) 
At ■ 



But since AW is of order v At, it then follows that 



AI4^^^/ 1 



At 



(27) 



so that dW/dt = oo as At ^ 0. 

Although the derivative of W{t) does not exist as a 
regular stochastic process, it is possible to give a math- 
ematical meaning to dW/dt as a generalized process (in 
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W(t) 50 



FIG. 7: Self-similarity of a Brownian motion path. In (a) we 
plot a path of a Brownian motion with 15000 time steps. The 
curve in (b) is a blow-up of the region delimited by a rectangle 
in (a), where we have rescaled the x axis by a factor 4 and 
the y axis by a factor 2. Note that the graphs in (a) and (b) 
"look the same," statistically speaking. This process can be 
repeated indefinitely. 



the sense of generalized functions or distributions). In 
this case, the derivative of the W{t) is called the white 
noise process ^(t): 



(28) 



I shall, of course, not attempt to give a rigorous definition 
of the white noise, and so the following intuitive argu- 
ment will suffice. Since according to H27|) the derivative 
dW 1 . , 

— ; — diverges as , a simple power-counting argument 

at ./dt 
suggests that integrals of the form 



m = / 9{mt')dt\ 



(29) 



should converge (in some sense); see below. 

In physics, the white noise £{t) is simply 'defined' as a 
'rapidly fluctuating function' [13 (in fact, a generalized 
stochastic process) that satisfies the following conditions 



m)) = 0, 
mm) = 5{t-t'). 



(30) 
(31) 



These two relations give the 'operational rules' from 
which quantities such as the mean and the variance of the 



integral I{t) in H29|l can be calculated. It is convenient, 
however, to have an alternative definition of stochastic 
integrals in terms of regular stochastic process. Such a 
construction was first given by the Japanese mathemati- 
cian Ito . 



C. Ito stochastic integrals 

Using (|28|l . let us first rewrite integral H29|) as an "in- 
tegral over the Wiener process" W{t): 



m = / git')dW{t'). 
Jo 



(32) 



The idea then is to define this integral as a kind of 
Riemann-Stieltjes integral. We thus take a partition 
{ti}i=Q of the interval [0, t] and consider the partial sums 



In = Y.9it^-i)^W{t,) 

i=l 



i=l 



g{t,-i)[W{t,)-W{t,.i)]. 

(33) 

The function g{t) above must satisfy certain appropriate 
conditions 0, the most important one being that g{t) be 
a non- anticipating function. This means, in particular, 
that the value g{ti-i) in (|33|) is independent of the 'next 
increment' /S.W{ti) of the Brownian motion. [For this 
reason, choosing to evaluate g{t) at the beginning of the 
interval Ai^ — ti — i^-i is a crucial point in the definition 
of the Ito stochastic integral. Another possible choice 
is to evaluate g{t) at the mid point t* = (ti-i + ti)/2, 
which leads to the Stratonovich integral 8]. In these 
notes I shall only consider Ito integrals.] 

Under the appropriate conditions on g{t), it is then 
possible to show that the partial sums /„ converge in the 
mean square sense. That is, there exists a process I{t) 
such that 



E {In-I{t)Y 



as n 



00. 



(34) 



Using the fact that g(t) is non-anticipating and that 
E [AW{t)] — 0, it follows immediately from the definition 
(j^ that /(t) has zero mean: 



E[I{t)] = E 



g{t')dW{t') 



= 0, 



(35) 



It is also possible to show that stochastic integrals obey 
the so-called isometry property: 



E 



{iit)y 



= E 



t ^ 2- 

g{t')dW{t') 



= / E[g^it')]dt'. 
Jo 

(36) 

We now see that the true meaning of conditions l|3U|) 
and H31|) is given by properties (|35|l and H3t)|) . for the 
particular case when g{t) is a deterministic function. 
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The Ito integral does not conform to the usual integra- 
tion rules from deterministic calculus. An example is the 
formula below 

^ WdW = ^W{t)' - ii, 

which is left as an exercise for the reader ltd inte- 
grals offer however a convenient way to define (and deal 
with) stochastic differential equations, as we will see next. 

D. Stochastic differential equations 

Physicists are quite familiar with differential equations 
involving stochastic terms, such as the Langevin equation 



where the constant fi represents the mean drift velocity. 
Integrating 1421) immediately yields the process 

X{t) ^ Ht + W{t), (43) 

whose pdf is 

p(..„._L.exp{(£^}^ ,44, 

Another important example of a (linear) SDE that can be 
solved explicitly is the geometric Brownian motion that 
will be discussed shortly. But before doing that, let us 
discuss a rather useful result known as ltd lemma or ltd 
formula. 

E. Ito formula 



dv 



(37) 



which describes_the motion of a Brownian particle in a 
viscous liquid JT^. Here 7 is the viscosity of the fluid 
and a is the 'amplitude' of the fluctuating force acting 
on the Brownian particle. (These parameters are usu- 
ally considered to be constant but in general they could 
be non-anticipating functions of time.) Equation H37|l 
does not however make much mathematical sense, since 
it evolves a quantity, namely, the derivative S^{t) of the 
Brownian motion, that does not even exist (except as a 
generalized process). Nevertheless, it is possible to put 
this equation on a firm mathematical basis by expressing 
it as a stochastic integral equation. First we rewrite H37fl 
as 



dv = — jvdt + adW, 
which upon integration yields 



(38) 



v{t) = u(0) - [ -fv{t')dt' + [ adW{t'). (39) 
Jo Jo 

This integral equation now makes perfectly good sense — 
in fact, its solution can be found explicitly [l9|. 

Let us now consider more general stochastic differential 
equations (SDE) of the form 



dX = a{X, t)dt + b{X, t)dW, 



(40) 



where a(x,t) and b{x,t) are known functions. Note that 
this 'differential equation' is actually a short-hand nota- 
tion for the following stochastic integral equation 

X{t)^X{0)+ [ a{X,t')dt' + [ b{XX)dW{t'). (41) 



Under certain condition on the functions a{x,t) and 
B{x, t), it is possible to show (see, e.g., @) that the SDE 
(|in|l has a unique solution X{t). 

Let us discuss another simple SDE, namely, the Brow- 
nian motion with drift: 



dX = fidt + adW, 



Consider the generic process X{t) described by the 
SDE l|in|) . and suppose that we have a new stochastic 
process Z defined by 

Z{t) = FiX{t),t), (45) 

for some given function F{x, t). We now wish to find the 
local dynamics followed by the Z{t), that is, the SDE 
whose solutions corresponds to the process Z(t) above. 
The answer is given by the Ito formula that we now pro- 
ceed to derive. 

First, consider the Taylor expansion of the function 

Fix,ty. 



-{dtf + -^7-7^dtdX 



2dtdx 



(46) 



a^idtf 



(47) 



2 dt^ 
Note, however, that 

{dXy = b'^dW'^ + 2abdtdW 
= b'^dt + 0{dt^/^), 

where we used the fact that dW^ = dt and dtdW = 
0{dt^^^). (Here we have momentarily omitted the ar- 
guments of the functions a and 6 for ease of notation.) 
Inserting H47() into 146() and retaining only terms up to 
order dt, we obtain 



dF = 



1 2 

'dt^2 'd^ 



dF 

dt + b -T^dX, 
ox 



(48) 



which is known as Ito formula. Upon using H40|) in the 
equation above, we obtain Ito formula in a more explicit 
fashion 



dF 



dF 
'dt 



a{X,t) 



dF 
dx 



1 d^F 



dt 



dF 

b{X,t)—dW, 



(49) 



What is noteworthy about this formula is the fact that 
the fluctuating part of the primary process X{t) con- 
tributes to the drift of the derived process Z{t) = F{t, X) 
through the term ^b'^{t,X)^^. We shall next use Ito 
(42) formula to solve explicitly a certain class of linear SDE's. 
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F. Geometric Brownian motion 

A stochastic process of great importance in Finance is 
the so-cahed geometric Brownian notion, which is defined 
as the solution to the foUowing SDE 



dS = ^iSdt + aSdW, 



(50) 



where fi and a are constants, subjected to a generic initial 
condition S{to) = Sq. Let us now perform the following 
change of variables Z = InS. Applying Ito formula 149|l 
with a = fiS, b = aS and F{S) = InS, it then follows 
that 



dZ 



1 



a^]dt + adW, 



(51) 



which upon integration yields 

Z{t) = Zo+(^^i-^a^^ {t-to) + a[Wit)-Wito)], (52) 

where Zq — In 5*0. Reverting to the variable S we obtain 
the explicit solution of the SDE jSUJl: 

S{t) = 5*0 exp I (^/i - ia^^ {t - <o) + <j[W{t) - W{to)] 

(53) 

From (|52|l we immediately see that Z{t) — Zq is dis- 
tributed according to A/" ((/i — jf^) t, (Jy/r^, where r = 
t—tQ. It then follows that the geometric Brownian motion 
with initial value S'(to) — Sq has the following log-normal 
distribution: 



where the time dependence has been omitted for clar- 
ity. We shall also often suppress the subscript from Vg{t) 
when there is no risk of confusion about to which port- 
folio we are referring. 

A portfolio is called self -financing if no money is taken 
from it for 'consumption' and no additional money is in- 
vested in it, so that any change in the portfolio value 
comes solely from changes in the asset prices. More pre- 
cisely, a portfolio x is self-financing if its dynamics is 
given by 



p(5,t; So, to) 



1 



V2^S 



exp ■ 



2 ^ 



(54) 

The geometric Brownian motion is the basic model for 
stock price dynamics in the Black-Scholes framework, to 
which we now turn. 



IV. THE STANDARD MODEL OF FINANCE 



dVsit) = x{t) ■ dS{t), t > 0. 



(56) 



The reason for this definition is that in the discrete-time 
case, i.e., t — tn, n — 0,1,2,..., the increase in wealth, 
^V{tn) = V{tn+i) — V{tn), of a self-finauciug portfolio 
over the time interval tn+i — t„ is given by 



AV{tn) = X{tn) ■ AS{tn), 



(57) 



where AS{tn) = S(tn+i) — S{tn). This means that over 
the time interval tn+i — tn the value of the portfolio varies 
only owing to the changes in the asset prices themselves, 
and then at time tn+i re-allocate the assets within the 
portfolio for the next time period. Equation (|56|l gen- 
eralizes this idea for the continuous-time limit. If fur- 
thermore we decide on the make up of the portfolio by 
looking only at the current prices and not on past times, 
i.e., if 

x{t) = x{t,S{t)), 

then the portfolio is said to be Markovian. Here we shall 
deal exclusively with Markovian portfolios. 

As we have seen already in Sec. Ill Dl an arbitrage rep- 
resents the possibility of making a riskless profit with no 
initial commitment of money. A more formal definition 
of arbitrage is as follows. 

Definition 4 An arbitrage is a portfolio whose value 
V(t) obeys the following conditions 

(i) y(o) = 

(a) V{t) > with probability 1 for all t > 
(Hi) V{T) > with positive probability for some T > 0. 



A. Portfolio dynamics and arbitrage 

Consider a financial market with only two assets: a 
risk-free bank account B and a stock S. In vector no- 
tation, we write S{t) = {B{t),S{t)) for the asset price 
vector at time t. A portfolio in this market consists of 
having an amount xq in the bank and owing xi stocks. 
The vector x{t) — {xQ{t),xi{t)) thus describes the time 
evolution of your portfolio in the {B, S) space. Note that 
Xi < means a short position on the ith asset, i.e., you 
'owe the market' \xi\ units of the ith asset. Let us denote 
by Vs{t) the money value of the portfolio x{t): 

Vs = X ■ S = xqB + xiS, (55) 



The meaning of the first condition is self-evident. The 
second condition says that there is no chance of losing 
money, while the third one states that there is a possi- 
bility that the portfolio will acquire a positive value at 
some time T. Thus, if you hold this portfolio until this 
arbitrage time there is a real chance that you will make 
a riskless profit out of nothing. [If P{V{T) > 0) = 1 we 
have a strong arbitrage opportunity, in which case we are 
sure to make a profit.] As we have already discussed in 
Sec. Ill Dl arbitrage opportunities are very rare and can 
last only for a very short time (typically, of the order of 
seconds or a few minutes at most). In fact, in the fa- 
mous Black-Scholes model that we will now discuss it is 
assumed that there is no arbitrage at all. 
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B. The Black-Scholes model for option pricing 

The two main assumptions of the Black-Scholes model 
are: 

(i) There are two assets in the market, a bank account 
B and a stock S, whose price dynamics are gov- 
erned by the following differential equations 



dB = rBdt, 

dS = ^iSdt + aSdW, 



(58) 
(59) 



where r is the risk-free interest rate, /i > is the 
stock mean rate of return, cr > is the volatil- 
ity, and W{t) is the standard Brownian motion or 
Wiener process. 

(ii) The market is free of arbitrage. 

Besides these two crucial hypothesis, there are addi- 
tional simplifying (technical) assumptions, such as: (iii) 
there is a liquid market for the underlying asset S as well 
as for the derivative one wishes to price, (iv) there are no 
transaction costs (i.e., no bid-ask spread), and (v) im- 
limited short selling is allowed for an unlimited period 
of time. It is implied by (|58|l that there is no interest- 
rate spread either, that is, money is borrowed and lent 
at the same rate r. Equation (|59|l also implies that the 
stock pays no dividend. [This last assumption can be 
relaxed to allow for dividend payments at a known (i.e., 
deterministic) rate; see, e.g., Q for details.] 

We shall next describe how derivatives can be 'ratio- 
nally' priced in the Black-Scholes model. We consider 
first a European call option for which a closed formula 
can be found. (More general European contingent claims 
will be briefly considered at the end of the Section.) Let 
us then denote by C{S, t; K, T) the present value of a 
European call option with strike price K and expiration 
date T on the underlying stock S. For ease of notation 
we shall drop the parameters K and T and simply write 
C{S, t). For later use, we note here that according to Ito 
formula (|49|l . with a — fiS and b — aS, the option price 
C obeys the following dynamics 



position on the option and a short position on A stocks. 
The value Il{t) of this portfolio is 

U{t) ^C{S,t)-AS. 

Since the portfolio is self-financing, it follows from H56|) 
that n obeys the following dynamics 



dIl = dC-A dS, 
which in view of (|59|l and H60() becomes 



(61) 



dn 



dC_ 

dt 



dS 



fiAS 



dt 



aS{^-A]dW. 



(62) 



We can now eliminate the risk [i.e., the stochastic term 
containing dW] from this portfolio by choosing 



dS 

Inserting this back into H62I) . wc then find 



dC 
'dt 



dt. 



(63) 



(64) 



Since we now have a risk-free (i.e., purely deterministic) 
portfolio, it must yield the same rate of return as the 
bank account, which means that 



dn = rUdt. 



(65) 



Comparing (|64|l with l|65() and using (|61|l and (|63() . we 

then obtain the Black-Scholes equation: 



dC 1 2q29^C 



rS—^rC^O, (66) 



which must be solved subjected to the following bound- 
ary condition 



C{S,T) = iascyi{S - K,Q). 



(67) 



dC 



dC 
'dt 



dC 1 



dt 



asf^dW. 



(60) 

In what follows, we will arrive at a partial differential 
equation, the so-called Black-Scholes equation (BSE), for 
the option price C{S, t). For pedagogical reasons, we will 
present two alternative derivations of the BSE using two 
distinct but related arguments: i) the A-hedging portfo- 
lio and ii) the replicating portfolio. 



The solution to the above boundary- value problem can be 
found explicitly (see below), but before going into that it 
is instructive to consider an alternative derivation of the 
BSE. [Note that the above derivation of the BSE remains 
valid also in the case that r, ^, and, a are deterministic 
functions of time, although a solution in closed form is 
no longer possible.] 



2. The replicating portfolio 



Ttie delta-hedging portfolio 



As in the binomial model of Sec. Ill El we consider the 
self-financing A-hedging portfolio, consisting of a long 



Here we will show that it is possible to form a portfolio 
on the {B,S) market that replicates the option C{S,t), 
and in the process of doing so we will arrive again at the 
BSE. Suppose then that there is indeed a self-financing 
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portfolio x{t) = {x{t),y{t)), whose value Z{t) equals the 
option price C{S,t) for all time t <T: 



Z EE xB + yS ^ C, 



(68) 



where we have omitted the time-dependence for brevity. 
Since the portfolio is self-financing it follows that 

dZ = xdB + yds = [rxB + fiyS)dt + aySdW. (69) 

But by assumption we have Z — C and so dZ = dC . 
Comparing 169() with (|60|l and equating the coefficients 
separately in both dW and dt, we obtain 



y ■ 



dC 

as' 



dc 1 , ,d^c 

Now from and {TUJl we get that 



1 

B 



C-S 



dC 

as 



(70) 
(71) 

(72) 



which inserted into H71|l yields again the BSE H66|) . as the 
reader can easily verify. 

We have thus proven, by direct construction, that the 
option C can be replicated in the (-B, S')-market by the 
portfolio {x,y), where x and y are given in l|72|l and H7()|l . 
respectively, with option price C being the solution of the 
BSE (with the corresponding boundary condition). [To 
complete the proof, we must also show that the initial 
price Co = C{S, 0) is the 'correct' one, in the sense that 
if the option price were Cq ^ Co, then there would be an 
arbitrage opportunity. In fact, if Cg > Cq an arbitrageur 
should short the option and invest in the replicating port- 
foho, whereas if Cq < Cq he should do the opposite.] 



The Black-Scholes formula 



After a somewhat tedious but straightforward algebra 
iSi, one obtains that in the new variables equation H66|) 
reads 



du d^u 
dr dx"^ ' 



(76) 



while the terminal condition H67f) becomes an initial con- 
dition 

u{x, 0) = uo{x) = max (e''^ - e"^, O) . (77) 

We now recall that the Green's function for the heat 
equation is 



C{x,x') = 



1 



' 4ttt 



SO that its generic solution for an arbitrary initial condi- 
tion uq{x) is given by 

/•OO 

u{x,t) = / uo{x')G{x,x')dx' 



1 



'Attt 



Uo 



(a;')e-("-"')'/^"dx'. (78) 



Inserting 1)77(1 into the integral above we obtain 
1 



u{t, x) 



V 4ttt Jo 



where 



/(a) 



/47rT 



(79) 



-(x-:.')V4r^^/^ (80) 



After completing the squares and performing some sim- 
plification, we find that 



Here we will solve equation (|66|l subjected to the 
boundary condition 16 7|) . Following the original work 
of Black and Scholes 12], the idea is to perform a change 
of variables so as to turn the BSE into the heat equation, 
which we know how to solve. Here we will not use the 
original transformation employed by these authors but a 
related one as shown below: 



T = 



T-t 



In 



where 



K 



1 /^2r \ ^ 1 /2r , 



(73) 



(74) 



(75) 



I{a) = 6'^^+'^ ^Nida), 



where 



da = 



X + 2aT 



(81) 



(82) 



and N(x) denotes the cumulative distribution function 
for a normal variable Af{0, 1): 



N{x) = 



1 



/27r 



'/'ds. 



(83) 



Inserting 1)81(1 into 1(79(1 and reverting back to the origi- 
nal dimensional variables, we obtain the famous Black- 
Scholes formula for the price of a European call option: 



C{S,t) = SN{di) - i^e-''(^-*)iV(d2), 



(84) 



14 



where 



di 



do 



(85) 
(86) 
(87) 



This formula is so often used in practice that it is aheady 
pre-defined in many software packages (e.g., Excel, Mat- 
lab, Maple, etc) as well as in most modern hand calcula- 
tors with financial functions. It should noted, however, 
that many people (academics and practitioners alike) be- 
lieve that the Black-Scholes model is too idealized to de- 
scribe real market situations; see Sees. V and VII for a 
brief discussion of possible extensions of the BS model. 



D. Completeness in the Black-Scholes model 

We have seen above that it is possible to replicate a 
European call option C{S, t) using an appropriate self- 
financing portfolio in the (-B, S) market. Looking back 
at the argument given in Sec. IIV B 21 we see that we 
never actually made use of the fact that the derivative 
in question was a call option — the nature of the deriva- 
tive appeared only through the boundary condition (|67fl . 
Thus, the derivation of the BSE presented there must 
hold for any contingent claim! 

To state this fact more precisely, let F{S, t) represent 
the price of an arbitrary European contingent claim with 
payoff F{S, T) = <&(S'), where $ is a known function. Re- 
tracing the steps outlined in Sec. lIVB^ we immediately 
conclude that the price ^(5*, t) will be the solution to the 
following boundary-value problem 



OF 
'dt 



1 2^2^^^ 



rS 



dF 



dS 
F{S,T) 



rF ^0 



(88) 
(89) 



Furthermore, if we repeat the arguments of preceding 
subsection and transform the Black-Scholes equation H88|) 
into the heat equation, we obtain that F{S, t) will be 
given by 



F{S,t) 



$(a;')e-(^-^')'/''^dx'. 



(90) 



where $(a;) denotes the payoff fmiction in terms of the 
dimensionless variable x; see (|73f) . Expressing this re- 
sult in terms of the original variables S and t yields a 
generalized Black-Scholes formula 



self-financing portfolio on the primary assets. Our 'proof 
of completeness' given above is, of course, valid only for 
the case of European contingent claims with a simple 
payoff function $(5'); it does not cover, for instance, 
path-dependent derivatives. It is possible however to give 
a formal proof that arbitrage-free models, such as the 
Black-Scholes model, are indeed complete; see Sec. IV CI 
Comparing the generalized Black-Scholes formula H91() 
with the pdf of the geometric Brownian motion given in 
(|54|l . we see that the former can be written in a conve- 
nient way as 



F{S,t)^e-^^'^-'^Efs[<^{ST% 



(92) 



F{S,t) 



-r{T-t) 



^27rcr2(T-t) 



$(5' 



(91) 



In summary, we have shown above that the Black- 
Scholes model is complete. A market is said to be com- 
plete if every contingent claim can be replicated with a 



where E^g[-] denotes expectation value with respect to 
the probability density of a geometric Brownian motion 
with fi = r, initial time final time T, and initial value 
S; see (|54|l . In other words, the present value of a con- 
tingent claim can be computed simply as its discounted 
expected value at maturity, under an appropriate proba- 
bility measure. This idea will become more clear after we 
discuss the notion of an equivalent martingale measure. 



V. EFFICIENT MARKETS: THE MARTINGALE 
APPROACH 

A. Martingales 

The concept of a martingale plays a important role 
in finance 20] . Unfortunately, a proper introduction to 
martingales req uires some knowledge of probability mea- 
sure theory |21| . Here however we will give a rather intu- 
itive discussion of martingales. For completeness we have 
listed in Appendix A some basic concepts from probabil- 
ity theory that would be required to make the following 
discussion more rigorous. 

We begin by recalling that a probability space is a triple 
{fl,J',P), where 

• is the space of elementary events or outcomes lo. 

• is a properly chosen family of subsets of SI, (i.e., 
a cr-algebra on SI). 

• P is a probability measure on T . 

In Finance, an outcome w is a 'market situation.' 
The family of subsets T specifies the class of events 
to which probabilities can be assigned. This is done 
through the concept of a cr-algebra, whose formal defi- 
nition is given in Appendix A. A probability measure P 
on T is simply a function P : — > [0, 1], satisfying a 
^•^e w 'obvious requirements': P(0) = 0, P(S7) = 1, and 
S'P{A^\JA2) = p\a^)+P{A2) if AinAa = 0. An element 
A oi T ^ A ^ T ^ is called a "measurable set" or "ob- 
servable event," meaning that it is possible to assign a 
"probability of occurrence," P(A) e [0, 1], for this event. 
Hence, T is said to be the set of 'observable events.' 
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Suppose now that we have a random function X : — > 
R. If to every subset of O of the form {lo : a < X{uj) < b} 
there corresponds an event A C then the function 
X is said to be measurable with respect to or simply 
!F-measurable. What this means is that it is possible to 
'measure' (i.e., assign a probability to) events of the form 
{a < X < b} through the obvious definition: P{{a < 
X < b}) = p{A). A J^-measurable function X is called a 
random variable. 

Let us next consider the notion of an "information 
flow." In a somewhat abstract way, we will represent 
the information available to an observer up to time t as 
a cr-algebra J-'t C J-^. In the context of Finance, the in- 
formation set J-'t would contain, for instance, the price 
history up to time t of all assets in the economy. It is 
natural to assume that J^s C J-t for s < t, since we expect 
that, as time goes on, we gain new information (and do 
not discard the old ones). Such a collection of ct- algebras 
represents an "information flow," or, more technically, a 
filtration. 

Definition 5 A filtration or information flow is a col- 
lection {J-t}t>o of a -algebras Tt T such that 

TsCTt, for 0<s<t. 

Suppose now we have a stochastic process Xt defined 
on (il,JF, P). (Think of Xt as being, say, the price of a 
given stock.) If the values of Xt can be completely de- 
termined from the information J^t, for alH > 0, then the 
process Xt is said to be adapted to the filtration {J-^t}t>a- 

Definition 6 The process Xt is adapted to the filtration 
{J't}t>o if Xt is Tt-measurable for all t>Q. 

A stochastic process Xt naturally generates an infor- 
mation flow, denoted by , which represents the "infor- 
mation contained in the trajectories of X{t) up to time 
t." A process X{t) is obviously adapted to its natural 
filtration . 

A last piece of mathematics is necessary to define a 
martingale, namely, the notion of conditional expecta- 
tion. This appears in connection with the crucial ques- 
tion of how the information, J-t^ , available at present time 
influences our knowledge about future values of X{t). If 
X(t) and Tto are not independent, then it is reasonable 
to expect that the information available up to the present 
time reduces the uncertainty about the future values of 
X{t). To reflect this gain of information is useful to in- 
troduce a new stochastic process 

Z{t) ^ E[Xt\Tt,i t>to, 

where the symbol E[Xt\!Fto] represents "the expected 
value of Xt, contingent on the information gathered up 
to time tg." 

A precise definition of conditional expectation is be- 
yond the scope of these notes. Here it will suffice to say 
that given a random variable 1" on a probability space 
(ri, P) and another a-algebra T' C it is possible 



to define a random variable Z = E\Y\T'], which rep- 
resents "the expected value of Y , given the information 
contained in JF'." The variable Z is a coarser version of 
the original variable Y, in the sense that we have used 
the information on J-' to reduce the uncertainty about Y . 
The following two properties of conditional expectations 
will be necessary later: 

E[Y\J^t] =Y, if y is J^t-measurable. (93) 
E[E[Y\J^t]]^E[Y]. (94) 

The first property above is somewhat obvious: if Y is 
JFj -measurable then Y and J^t 'contain the same infor- 
mation,' hence taking expectation of Y conditional to J^t 
does not reduce the uncertainty about Y. The second 
property is the so-called law of iterated expectations, and 
basically represents the law of total probability. 

After these mathematical preliminaries, we are now in 
a position to define martingales. 

Definition 7 A stochastic process Mt is called a martin- 
gale with respect to the filtration {J^t}t>o if 

(i) Mt is adapted to the filtration {J-t}t>o 

(ii) E[\Mt\] < oo for all t > 

(ill) E[Mt\Tt,] = Mt, for all t > to 

Condition (i) simply says that Mt can be determined 
from the information available up to time t, whereas con- 
dition (ii) is a technicality. The defining property of a 
martingale is therefore condition (iii), which is usually 
referred to as the martingale condition. It says that the 
best prediction of future values of the process M{t), con- 
tingent on the information available at the present time, 
is the current value Mt^ . 

Because of property (iii), a martingale is usually de- 
scribed as a "fair game." To see why this is so, suppose 
that Mn represents the fortune of a gambler at time n. 
(For convenience let us assume here that time is discrete.) 
The difference /i„ = Af„ — M„_i is then the amount the 
gambler wins on the nth play (a negative win is of course 
a loss). Now let us compute the gambler's expected gain 
on the (n -f l)th play, given the information up to time 
n: 

E[hr^+l\J^n] = E[M,,+ i - MnlJ'n] 

= £;[M„+i|.F„]-S[M„|.F„] 

= Mn-Mn 

= 0, (95) 

where in the third equality we used the martingale prop- 
erty and rule H93|l . We thus have that at each new play 
of the game the expected gain is null, and in this sense 
it is a "fair" game. 

A Brownian motion W{t) is a martingale with respect 
to its natural filtration. To show this, we only need to 
verify the martingale condition (since the other two con- 
ditions are trivially fulfilled): 

E[W{t)\Tt,] = E[W{t)-W{to) + W{to)\Tt,] 
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= E[W(t)-W{to)\Tu,] 
= + Wito) 



E[W{to)\Tu,] 
(96) 



In the third equahty above we used the fact that the in- 
crements W{t) — W{to) are independent of J^to and have 
zero mean, together with property (|93|l . It is also possi- 
ble to show that Ito stochastic integrals are martingales. 
Indeed, the theory of stochastic integration is intimately 
connected with martingale theory [5j . 

Another important property of martingales is that 
their expected value remains constant in time: 

E[Mo] ^ E[E[Mt\To]] = E[Mt], 

where in first equality we used the martingale property, 
while in the second equality property 194(1 was used. 
Thus, a necessary (but not sufficient) condition for a pro- 
cess to be a martingale is that it have no drift. Therefore, 
a diffusion process of the form 14U|) is not a martingales 
unless the drift term vanishes. For this reason, the geo- 
metric Brownian motion H5U|) is not a martingale. It is 
possible, however, to introduce a new probability mea- 
sure Q, with respect to which the geometric Brownian 
motion becomes a standard Brownian motion and hence 
a martingale, as discussed next. 



B. Equivalent martingale measures 

Recall that a probability measure P on a measurable 
space (ri, is a function P : T ^ [0,1] that assigns to 
every event A C J- a real number P{A) € [0, 1]. Suppose 
now we have another probability measure Q defined on 
the same space {il,J-). We say that the probability mea- 
sures P and Q are equivalent if the following condition 
is satisfied: 



Q{A) = ^ P{A) = 0, for aU AeT. 



(97) 



To get a better grasp on the meaning of the condition 
above, consider a random variable X [on P)]. If Q 

is a probability measure equivalent to P, then condition 
(|97|l implies that there exists a function p{X) such that 
expected values w.r.t Q are calculated in the following 
way 



EQ[g{X)]=Ep[p{X)g{X)], 



(98) 



where g{x) is an arbitrary function. Alternatively, we 
may write (|98|) in terms of probability densities: 



fQ{x) = p{x)fp{x), 



(99) 



where fp{x) denotes the probability density of X 
w.r.t the measure P and fgix) is the density w.r.t Q. 
(The function p{x) is called the Radon-Nikodym deriva- 
tive of measure Q with respect to measure P ^Y\.) 
Consider now the Brownian motion with drift 



W{t) ^at + W{t), 0<t<T 



(100) 



where a is some constant. (The finite-horizon condition 
t < T is a technicality that is not relevant for our pur- 
poses.) As already noted, W{t) is not a martingale (since 
its expected value is not constant). However, there is an 
equivalent probability measure Q, with respect to which 
the process W{t) becomes the standard Brownian mo- 
tion (and hence a martingale). This result is known as 
Girsanov theorem. 

Theorem 1 (Girsanov theorem) The process W{t) 
given in 1^1 UUj) is a standard Brownian motion with re- 
spect to the probability measure Q defined by 



fQ{x,t) ^ Mt{i)fp{i,t), 
where Mt is the process 
1 



(101) 



1 



Mt = exp i^-aWt - tj = cxp |-aM^t -I- -a^t 

(102) 

Proof. For a formal proof see, e.g., 0. Here we shall only 
sketch a proof of the fact that the process W{t) is indeed 
distributed according to A/'(0, Vt), as the standard Brow- 
nian motion. First recall from (|44|) that the probability 
density fp{x, t) of W{t) under the original measure P is 



/P(.,t) = -=exp| ^ 



Now according to (|101|l and (|102|l we have that 



/Q(£,t) = e-''^+^'*/p(i,0- 



Inserting H103|) into (|104|l then yields 



V2t 



(103) 



(104) 



(105) 



which is precisely the pdf for the standard Brownian mo- 
tion. Q.E.D. 

One of the main applications of change of measures is 
to eliminate the drift in stochastic differential equations, 
so that with respect to the new measure Q the process is 
a martingale. The measure Q is then called an equivalent 
martingale measure. Constructing the equivalent martin- 
gale measure for an arbitrary SDE of the form H41|l is a 
rather complicated procedure |^. One important excep- 
tion are linear SDE's where the Girsanov theorem gives 
the measure transformation in an explicit form, as shown 
below. 

Consider the geometric Brownian motion discussed in 
Sec. IIIIFI For technical reasons 8], let us restrict our- 
selves to its finite-horizon version: 



dS = pSdt + aSdW, t < T. 



(106) 



where p and a are positive constants. This equation can 
then be rewritten as 



dS = aS (^dt + dW] = aSdW, 



(107) 
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where 



Wt^{li/a)t + Wu t<T. 



(108) 



eQ 



Now, according to Girsanov theorem, Wt is a standard 
Brownian motion with respect to the measure Q given 
in (|101|l with a — /i/cr, and since the SDE (|107|l has no 
drift, its solution St is a martingale w.r.t. the measure 



The 'efficiency symmetry' and the no-arbitrage 
condition as its 'conservation law' 



The notion of an efficient market plays a crucial role 
in Finance. Roughly speaking, in an efficient market all 
relevant information is already reflected in the prices ■ 
This means, in particular, that past prices give no ad- 
ditional information that is not already contained in the 
current price. In an efhcient market prices thus adjust 
immediately to the arrival of new information. Since the 
content of future information and their effect on prices 
are unknown, it should be impossible to make definite 
predictions about future price based on the informa- 
tion available today. Thus, the best prediction for the 
expected future price (discounted to the present time) 
should be today's price. In other words, in a efficient 
market 'discounted prices' should be a martingale. This 
intuitive definition of efficiency is formalized below. [For 
technical reasons the results of this section will be re- 
stricted to the discrete-time case. Their extension to the 
continuous-time framework, although possible to a large 
extent, is more complicated and will not be attempted in 
the present notes; see, e.g., 0.] 

Consider a market formed by two assets {B,S), where 
B is our usual risk-free free asset and S* is a risky as- 
set. We suppose that S{t) follows a stochastic process on 
a probability space {Q,J^,P) endowed with a filtration 
{^t}t>o- 

Definition 8 Suppose the market {B, S) operates at dis- 
crete time tn, n = 1,2, .... This market is said to be effi- 
cient if there exists (at least one) probability measure Q, 



equivalent to P , such that the 'discounted price 



,s(ti 

B{t) 



a martingale with respect to the measure Q, that is, 



eQ 



s{t) 

B{t) 



to 



Sjto) 
B{to)' 



for to <t. 



(109) 



S{t + T) 
B{t + T) 



Bity 



(110) 



for any T > 0. 

In Physics, symmetry principles are intimately con- 
nected with conservation laws. (Recall, for instance, that 
the invariance of Newton's law under time translation im- 
plies conservation of energy.) It is thus only natural to 
ask whether the 'efficiency symmetry' above also leads 
to a 'conservation law.' Perhaps not surprisingly, this 
is indeed the stated in the following theorem, 

which is sometimes referred to as the First Fundamental 
Theorem of asset pricing. 

Theorem 2 Suppose the market {B, S) operates at dis- 
crete time tn, n = 1,2,.... Then this market is efficient 
if and only if it is arbitrage-free. 

(See 13 for a proof.) 

Recall that absence of arbitrage means that any self- 
financing portfolio with zero initial value, and with no 
chance of becoming negative, will remain zero- valued for 
all subsequent times. More precisely, the no-arbitrage 
condition says that 



V^(0) = and V{t) > a.s. 



V{t) = a.s., (Ill) 



where E'^ means expected value w.r.t. the measure 



where a.s. means almost surely, i.e., with probability 1. 
The absence of arbitrage can thus be interpreted as a 
kind of 'conservation law' for the "vacuum state" of the 
market: if you start in a state with zero initial money 
and do not take any risks, then you remain at this state 
for all times. I thus find it quite interesting that the so- 
called "no-free-lunch" condition can actually be seen as 
the conservation law associated with the efficiency sym- 
metry. 

Another very important result is the Second Funda- 
mental Theorem of asset pricing, linking the completeness 
of a market to the uniqueness of its equivalent martingale 
measure. 

Theorem 3 An arbitrage- free {B, S)-market is complete 
if and only if the equivalent martingale measure Q is 
unique. 

(See 13 for a proof.) 

We already know that the Black-Scholes model is com- 
plete. Below we will calculate its equivalent martingale 
measure explicitly, and it will become clear from the con- 
struction that it is indeed unique. 



The requirement of efficiency, as defined above, 
is somewhat reminiscent of a symmetry principle in 
Physics. Indeed, we can recast definition (|109|l by saying 
that in a efficient market there exists a 'special measure 
Q' with respect to which discounted prices are invariant 
under a sort of 'time translation,' in the following sense: 



D. Pricing derivatives with the equivalent 
martingale measure 

The notion of an equivalent martingale measure can be 
used to price derivatives in a rather direct way, without 
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having to solve a PDE. The idea is that in an efficient 
economy all financial assets are martingales with respect 
to the equivalent martingale measure Q. More precisely, 
if F{S^ t) is a contingent claim with maturity T and payoff 
function $(S'(T)) then from H109I) we have 



F{S,t) 
Bit) 



K 



t,s 



^{S[T)) 



B{T) 



(112) 



where the subscripts S denote that the expected value 
is taken at present time t and with current value S, i.e., 
conditional to the information available at time t. It is 
not hard to convince oneself that if the derivative price F 
were not given by H11U|) , then there would be an arbitrage 
opportunity. 

In the Black-Scholes model, the risk-free asset is a bank 
account with fixed interest rate r, i.e., B{t) = e''*, so that 
(I112|l becomes 

F{S,t)^e-^^'^-'^Ef,MS{T))], (113) 



or 



F{S,t) 



-r{T-t) 



<^{S')fQ{S\T-S,t)dS', (114) 



where fqiS, t; 5o, to) denotes the probability density, un- 
der the equivalent martingale measure Q, of the process 
S{t) with initial value S{to) — Sq. All that remains to 
be done now is to find the equivalent martingale measure 
Q for the Black-Scholes model. To do this, consider the 
process 



Zit) 



B{t) 



(115) 



We then have 



dZ = 



re -""^Sdt + e-^'^dS 
{fi-r)Zdt + (jZdW 
aZdW, 



where 



W{t) = [{^l-r)/a]t + W{t). 



(116) 



Now recall that in the Black-Scholes model the stock 
price follows a geometric Brownian motion 

dS = ^iSdt + aSdW, (117) 

which in terms of the process W{t) given in (|116|l reads 



dS = rSdt + aSdW. 



(118) 



From Girsanov theorem we know that there is a equiv- 
alent martingale measure Q that turns W{t) into a Brow- 
nian motion. Equation (|118|l then shows that w.r.t the 
measure Q the price S{t) follows a geometric Brownian 
motion with mean rate of return equal to r. The proba- 
bility density fQ{S', T; 5*, t) can now be obtained directly 



from (|54() , by simply setting /i 
gets 



fQ{S',T; S,t) 



1 



exp ■ 



r and 5*0 = 5*. One then 



2ah 



(119) 



where t = T — t. Inserting H119|l into (|114|l we obtain 



F{t,S) = 



,-r{T-t) 



$(5') 



,[ln(S7S)-(r-i^^)(T-t)]'^ 

S' 



(120) 

which is precisely the expression obtained for the generic 
solution of the Black-Scholes equation given in (|92|l . In 
the case of a European call option we have $(5") = 
max(S" — K,0), which inserted into H12U|) yields, after 
some algebra, the Black-Scholes formula H84|l . 

It is interesting to notice that under the equivalent 
martingale measure Q the stock price in the Black- 
Scholes model follows a geometric Brownian motion with 
the mean rate of return equal to the risk-free interest 
rate r; see H118|l . It is as if all investors were risk neu- 
tral, in the sense they would be willing to invest on a 
risky stock even though its expected return is just what 
a risk-free bank account would yield. For this reason, 
the pricing method based on the equivalent martingale 
measure is commonly referred to as risk neutral valua- 
tion. Of course, actual investors are not risk neutral. 
However, in an efficient market there is a 'special refer- 
ence frame' where investors can be treated as if they were 
indeed insensitive to risk. 



VI. BEYOND THE STANDARD MODEL OF 
FINANCE I: NON-GAUSSIAN DISTRIBUTIONS 

We have seen above that the Black-Scholes model is 
an elegant and powerful theoretical construct: it is com- 
plete, efficient, arbitrage-free, and Gaussian (in the sense 
that the stock returns are normally distributed). It is 
thus important to ask whether real markets actually fit 
into this nice framework. 

There are two main ways in which real markets may 
deviate from the standard Black-Scholes model: i) the 
returns may not be normally distributed or ii) there may 
exist long-memory effects on the time series. In this Sec- 
tion we will discuss the possibility that asset prices may 
follow a non-Gaussian stable Levy process, while in the 
next section we investigate whether financial data might 
exhibit long-term memory. 

Mandelbrot |2^ in 1963 was perhaps the first person 
to challenge the paradigm that returns are normally dis- 
tributed. He analyzed cotton prices on various exchanges 
in the United States and found evidences that their dis- 
tribution of returns decays as a power law and hence 
much slower than a Gaussian. An example of this 'fat 
tail' behavior can be seen in Fig. O where I plot the dis- 
tribution for the returns of the Ibovespa index. In this 
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FIG. 8: Distribution of the daily Ibovespa returns (open cir- 
cles). The solid line correspond to a Gaussian distribution 
with the same variance as that of the empirical distribution. 



figure, it is also shown a Gaussian distribution with the 
variance of the data, which appears as a parabola in the 
linear-log scale of the graph. One clearly sees that the 
empirical distribution does indeed have 'fatter tails' when 
compared with the Gaussian distribution. 

The Gaussian distribution is special for two main rea- 
sons. The first one is the Central Limit Theorem |23| 
that states that the sum of infinitely many independent 
random variables (with finite variance) will converge to 
a Gaussian variable. The second one is the fact that it 
is a stable distribution, in the sense that the sum of two 
independent Gaussian random variables is also a Gaus- 
sian variable. It is thus natural to ask whether there are 
other stable distributions. The French mathematician 
Paul Levy showed that there is indeed a whole family 
of stable distributions of which the Gaussian is but one 
particular case. In what follows, I will first introduce 
the so-called Levy stable distributions and then briefly 
discuss their possible applications to financial data. 

A. The stable Levy distributions 

Let X be a random variable with probability density 
function (pdf) given by p{x). We recall that the char- 
acteristic function ip{z) of a random variable X is the 
Fourier transform of its pdf p{x): 

/CO 
p{x)e''"'dx. (121) 
-oo 

Now let Xi and X2 be two independent random vari- 
ables with pdf's Pi{xi) and P2{x2)- Since Xi and X2 are 
independent, the pdf of the sum X = Xi + X2 is the 
convolution of the original pdf's: 

/oo 
Pi{s)p2{x - s)ds. 
'OO 

Let now ip{z, 2) denote the characteristic function of X. 
In view of the convolution theorem which states that the 



Fourier transform of the convolution is the product of the 
Fourier transforms, it then follows that 

^iz,2)^^i{z)^2{z). (122) 

Suppose furthermore that Xi and X2 are identically dis- 
tributed, that is, 

Xi = X2, (123) 

or alternatively 

^i{z) = ^2{z)^ip{z). (124) 

(Recall that the symbol ^ denotes equality in the dis- 
tribution sense.) From (|122|l and 1)124(1 it then follows 
that 

(^(z,2) = [^(z)]2. 

In general, if X — X^il^i where the X^'s are inde- 
pendent and identically distributed (i.i.d.) random vari- 
ables, then 

^{z,N)^Mz)f, (125) 

from which the pdf of X can be obtained by calculating 
the inverse Fourier transform. Note that the pdf of the 
sum of N i.i.d. random variables will in general be quite 
different from the pdf of the individual variables. There 
is however a special class of distribution, the stable dis- 
tributions, for which the pdf of the sum has the same 
functional form of the individual pdf's. 

Definition 9 A probability distribution p{x) is stable if 
for each N > 2 there exist numbers ajv > and bjy, 
such that, if Xi, Xn are i.i.d. random variables with 
distribution p{x) , then 

Xi + ... + Xn ^anXi + bN. (126) 

In other words, a distribution is stable if its form is 
invariant under addition, up to a rescaling of the variable 
by a translation and a dilation. More precisely, if p(x, N) 
denotes the probability density of X = X]i=i where 
the XiS are i.i.d. variables with a stable distribution 
p{x), then H12t)|) implies that 

p{x,N)^^p(^—^). (127) 

Stability of probability distributions has a nice interpre- 
tation from an economic standpoint. To see this, sup- 
pose that the variables Xi represent daily increments of 
a given financial asset. Stability then means preserva- 
tion of the distribution under time aggregation, which is 
a rather natural property to expect from financial data. 

As already mentioned, the Gaussian distribution is sta- 
ble. To see this, recall that the Fourier transform of a 
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Gaussian is a Gaussian and that the product of Gaus- 
sians is again a Gaussian, so that from (|125|) it follows 
that the characteristic function ip{z, N) will indeed be 
that of a Gaussian variable. More precisely, we have that 
the characteristic function of a normal variable M{0, a) 
is (p{z) — e~^°' /^^^ , which inserted into H125|) immedi- 
ately yields ip{z,N) = e~(^'^'/2)z" _ ^hus, the sum of N 
i.i.d. normal variables is normally distributed with stan- 
dard deviation y/Na, that is, X = N{Q,VNu), which in 
turn implies that 



p{x,N) = -)=P ( —7= 



(128) 



Thus, in the case of the Gaussian distribution we have 
ajv = and 5jv — 0- 

The class of stable distributions is rather small and was 
completely determined by the mathematicians P. Levy 
and A. Ya. Khintchine in the 1920's. Here we shall re- 
strict ourselves to the subclass of symmetric distribution. 
In this case, the characteristic function is given by 



'a\z\^ 



(129) 



where < a < 2 and a > 0. The parameter a is called 
the stability exponent and a is a scale factor. Taking the 
inverse Fourier transform of ifaiz) we obtain the corre- 
sponding pdf Paix): 



Ztt 



1 

dz ~ — I e 

TI" ./o 



cos{zx)dz . 
(130) 



Unfortunately, however, only for two particular values of 
a can the integral above be explicitly calculated: 

• a = 1 (Lorentzian or Cauchy distribution): 

2a 1 

p[x) 



IT X^ 



4a2 



a = 2 (Gaussian distribution): 



We thus see that the Levy distribution with a < 2 has 
the interesting property that it shows scaling behavior 
for large x, i.e., p{x) decays as a power-law. 

Another important scaling relation for the Levy dis- 
tribution can be obtained, as follows. First note that 
for symmetric stable distribution we necessarily have 
— 0. Now using l|125|l and (|129|l . we easily find that 
the dilation factor a at in (|127|l is 



OAT 



so that (|127(l becomes 



Paix,N) 

which implies that 



Pa (N^/'^x) 



p{0,N) 



ATI/"' 



(133) 



(134) 



(135) 



One can then use this scaling relation to estimate the 
index a of the Levy distribution: in a log-log plot of 
p(0, N) against N, the slope of a linear fit gives precisely 
1/a; see Sec. IVI CJl below. 

The power-law decay of Levy distributions implies, of 
course, the absence of a characteristic scale. The down- 
side of this is that all Levy distributions have infinite 
variance! In fact, all moments of order higher than 1 
are infinite, since diverges for n > a, as can be 

readily shown from H131|) . Processes with infinite vari- 
ance are not physically plausible, so several prescriptions 
to truncate the Levy distribution at some large scale has 
been proposed, as discussed next. 



B. Truncated Levy distributions 

To circumvent the problem of infinite variance in Levy 
distributions, several truncation prescriptions have been 
proposed in the literature. In a general they can be writ- 
ten as 



' Sira 

Note also that Levy distributions are not defined for 
a > 2, because in this case the function obtained from 
(|130|l is not everywhere positive. 

Although, for arbitrary a the pdf Pa{x) cannot be 
found in closed form, its asymptotic behavior for large 
X can be easily calculated from (|130|l . Here one finds 
that 



where 



Paix) 



|1 + Q ■ 



00, 



- r(l-l-a) sm— . 

TT 2 



p{x) ^Pa{x)<^{x), 



(136) 



where $(x) is a cut-off function to be chosen in such way 
that the variance of the truncated distribution is finite. 
For example, two possible choices that have been used to 
model the distributions of financial asset prices are given 
below 

• Abruptly truncated Levy distribution (ATLD): 



(131) ^x)^eixc-\x\), 

where Q{x) is the Heaviside function and Xc is some cut- 
off length scale. 

(132) 

• Exponentially truncated Levy distribution (ETLD): 
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where A > and A is a normalization factor. 

Other variants of truncated Levy distributions that 
have also been considered in the literature are the grad- 
ually truncated Levy distribution 1231 and the exponen- 
tially damped Levy distributions |24| . 

Since a truncated Levy distribution has finite variance, 
then by the central limit theorem the distribution of the 
sum X = Xi + ... + Xn of N i.i.d variables with such 
a distribution will converge to a Gaussian distribution 
for large N. However, this convergence is usually very 
slow — for financial data it is typically of the order of tens 
of days; see below. For shorter time scales, non-Gaussian 
behavior may thus be of practical relevance. 

C. Levy distributions in Finance 



for N < 20 days, whereas for larger N a Gaussian-like 
behavior (i.e., a = 0.5) was observed. 

Many other applications of Levy processes in Finance 
have been discussed in the literature |2^. For example, 
a model for option pricing has recently been considered 
where the price of the underlying asset is assumed to fol- 
low a truncated Levy process |27[. More recently, there 
have accumulated evidences |28l l29l |^] that in certain 
cases financial data may be better described by expo- 
nential distributions, rather than by Levy or Gaussian 
distributions. 



VII. BEYOND THE STANDARD MODEL OF 
FINANCE II: LONG-RANGE CORRELATIONS 

In this section, we discuss the possibility that asset 
prices might exhibit long-range correlations and thus may 
need to be described in terms of long-memory processes, 
such as the fractional Brownian motion. 



Levy distribution have been used, for example, by 
Mantegna & Stanley |^ to model the distribution of 
changes in the stock index S&P500 of the American 
Stock Exchange. They analyzed high-frequency data 
(one-minute quotes of the S&P500) over the period from 
January 1984 to December 1989. From the original time 
series Y{t) of the index values, they first generated time 
series corresponding to index changes during intervals of 
N minutes: 

ZN{t) = Y{t + N)~Y{t). (137) 

They then computed the empirical pdf p(z, N) and an- 
alyzed the scaling of p(0, N) with N. In a log-log plot 
p{0,N) showed a linear behavior, as predicted by (|135|l . 
with a slope corresponding to a = 1.4 for 30 < < 1000 
minutes For N > 10^ the slope of p(0, N) approaches 
—0.5, indicating convergence to a Gaussian behavior. In 
the Levy regime (i.e., small N), however, the tail of their 
empirical pdf decays slower than Gaussian but faster 
than a pure Levy distribution with the exponent a found 
from above scaling argument. These facts thus suggest 
that a truncated Levy distribution would perhaps be 
more appropriate for modeling the actual distribution. 
Indeed, Bouchaud and Potters |0 found that the prob- 
ability of 15-minute changes of the S&P500 index is well 
described by a ETLD with a = 1.5. 

The ETLD has also been applied by Miranda & Riera 
[25| to study the daily returns of Ibovespa index of the 
Sao Paulo Stock Exchange in the period 1986-2000. From 
the daily closing values Y{t) of the Ibovespa, they first 
calculated the time series for the returns in intervals of 
N days 

rN{t) = \ogY{t + N)-logY{t), (138) 

and then computed the corresponding pdf's for p{r,N). 
From the scaling of p(0,iV) they found a ~ 1.6 — 1.7 



A. Fractional Brownian motion 

The fractional Brownian motion (FBM) is a Gaussian 
process {Wif(t),i > 0} with zero mean and stationary 
increments, whose variance and covariance are given by 

E[Wl(t)]=f", (139) 

E[WH{s)WH{t)] = \ {s^" + f^" -\t- sp^)(J40) 

where Q < H < 1. The FBM Wnit) is a self-similar 
process, in the sense that 

WH{at)=a"WH{t), (141) 

for all a > 0. A sample path of a FBM is therefore a 
fractal curve with fractal dimension D = 1/H. The pa- 
rameter H is called the self-similarity exponent or the 
Hurst exponent. For H = 1/2 the process Wnit) corre- 
sponds to the usual Brownian motion, in which case the 
increments Xt = Wnit + 1) — Wnit) are statistically in- 
dependent, corresponding to white noise. On the other 
hand, for H ^ 1/2 the increments Xt, known as frac- 
tional white noise, display long-range correlation in the 
sense that 

E[Xt+hXt] ~ 2H{2H - l)h'^"-^ for h-^oo, (142) 

as one can easily verify from H139() and l|140|l . Thus, if 
1/2 < < 1 the increments of the FBM are positively 
correlated and we say that the process Wnit) exhibits 
persistence. Likewise, for < < 1/2 the increments 
are negatively correlated and the FBM is said to show 
antipersistence. Sample FBM paths with H = 0.2,0.5, 
and 0.8 are shown in Fig. 

Several estimators for the exponent H have been dis- 
cussed in the literature; see, e.g., Ref. for a com- 
parison among some of them. One general methodology 
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where S is the data standard deviation 



FIG. 9: Sample paths of fractional Brownian motion. 



consists in estimating how the 'amount of fluctuation' 
within a time window of size r scales with t. Specific 
methods, such as the Hurst rescaled range (R/S) analy- 
sis [s^ or the Detrendend Fluctuation Analysis j33j _34jJ : 
differ basically on the choice of the fluctuation measure. 
Here I shall discuss only the DFA that has proven to 
be a more reliable estimator for H than the Hurst R/S 
analysis (ssf. 



B. Detrended fluctuation analysis 

Suppose we have a time series r(t), t = 1,...,T, cor- 
responding to, say, daily returns of a financial asset. To 
implement the DFA, we first integrate the original time 
series r{t) to obtain the cumulative time series X{t): 



S 



\ 



1 ' 



n - r) 



(147) 



The Hurst exponent H is then obtained from the scaling 
behavior of F{t): 



F{t)=Ct", 



(148) 



where C is a constant independent of the time lag r. 

In a double-logarithmic plot the relationship 1] 148(1 
yields a straight line whose slope is precisely the expo- 
nent 7J, and so a linear regression of the empirical F{t) 
will immediately give H . One practical problem with this 
method, however, is that the values obtained for H are 
somewhat dependent on the choice of the interval within 
which to perform the linear fit [s^ [3^1 . It is possible to 
avoid part of this difficulty by relying on the fact that for 
the fractional Brownian motion^ the fluctuation function 
F(t) can be computed exactly [Slj : 



where 



Ch 



Fh{t) = Cht 



1 



H 



2H +1 H + 2 H + 1 



1/2 



(149) 



(150) 



In H149|) we have added a subscript H to the function 
F to denote explicitly that it refers to Wnit). Equation 
(|149|l with 1)150(1 now gives a one-parameter estimator 
for the exponent H: one has simply to adjust H so as to 
obtain the best agreement between the theoretical curve 
predicted by Fh{t) and the empirical data for F{t). 



X{t)^Y.^T{t')-r), t = l,...,T, 



t'=i 



where 



(143) 



(144) 



Next we break up X{t) into N non-overlapping time in- 
tervals, /„, of equal size r, where n = 0, l,...,iV — 1 and 
iV corresponds to the integer part of T/r. We then in- 
troduce the local trend function IV (0 defined by 



Yr{t)^an+bnt for t e In, 



(145) 



where the coefficients a„ and 6„ represent the least- 
square linear fit of X{t) in the interval /„. Finally, we 
coinpute the rescaled fluctuation function F{t) defined 
as Ha 



\ 



Nt 

-j2[x{t)-YAt)r 

TLT ^ — ' 



(146) 



C. Fractional Brownian motion in Finance 

The idea of using the FMB for modeling asset price 
dynamics dates back to the work of Mandelbrot & van 
Ness ■ Since then, the Hurst exponent has been cal- 
culated (using different estimators) for many financial 
time series, such as stock prices, stock indexes and cur- 
rency exchange rates [s^ |3^ [s^, E3j ^ • In many cases 
[ssj an exponent H > 1/2 has been found, indicating 
the existence of long-range correlation (persistence) in 
the data. It is to be noted, however, that the values of 
H computed using the traditional R/S-analysis, such as 
those quoted in , should be viewed with some caution, 
for this method has been shown j35| to overestimate the 
value of H. In this sense, the DFA appears to give a more 
reliable estimates for H. 

An example of the DFA applied to the returns of the 
Ibovespa stock index is shown in Fig. ^| (upper curve). 
In this figure the upper straight line corresponds to the 
theoretical curve Fh{t) given in ((149(1 with H = 0.6, and 
one sees an excellent agreement with the empirical data 
up to T ~ 130 days. The fact that H > 0.5 thus indicates 
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FIG. 10: Fluctuation function F(t) as a function of r for 
the returns of the Ibovespa index (upper curve) and for the 
shufhed data (lower curve). The upper (lower) straight line 
gives the theoretical curve Fh{t) for H — 0.6 (H — 1/2). 



persistence in the Ibovespa returns. For r > 130 the data 
deviate from the initial scaling behavior and cross over 
to a regime with a slope closer to 1/2, meaning that the 
Ibovespa looses its 'memory' after a period of about 6 
months. Also shown in Fig. ^| is the corresponding F{t) 
calculated for the shuffled Ibovespa returns. In this case 
we obtain an almost perfect scaling with H = 1/2, as 
expected, since the shuffling procedure tends to destroys 
any previously existing correlation. 

As already mentioned, the Hurst exponent has been 
calculated for many financial time series. In the case of 
stock indexes, the following interesting picture appears 
to be emerging from recent studies |35ll40ll4]| : large and 
more developed markets, such as the New York and the 
London Stock Exchanges, usually have H equal to (or 
slightly less than) 1/2, whereas less developed markets 
show a tendency to have H > 1/2. In other words, large 
markets seem indeed to be 'efflcient' in the sense that 
H ~ 1/2, whereas less developed markets tend to ex- 
hibit long-range correlation. A possible interpretation for 
this finding is that smaller markets are conceivably more 
prone to 'correlated fluctuations' and perhaps more sus- 
ceptible to being pushed around by aggressive investors, 
which may explain in part a Hurst exponent greater than 
1/2. 

It should also be pointed out that a considerable time- 
variability of the exponent H for stock indexes has been 
found lla, ^3 1 indicating that the data in such cases 
cannot be modeled in terms of stationary stochastic pro- 
cesses. (A time- varying H has also been observed in other 
financial data, such as, currency exchange rate |39||.) In 
such cases, the stationarity assumption is only a rather 
crude approximation |35j . Furthermore, the time de- 
pendence of the Hurst exponent is an indication that 
the underlying process might be multifractal rather than 
monofractal; see below. 

A time-varying Hurst exponent has been observed for 
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t 

FIG. 11: The Hurst exponent H for the Ibovespa as a function 
of time. Here H was computed in two-year time intervals, 
with the variable t denoting the origin of each such interval. 



the Brazilian stock market. This case is of particular in- 
terest because in the recent past Brazil was plagued by 
runaway inflation and endured several ill-fated economic 
plans designed to control it. To analyze the effect of 
inflation and the economic plans, Costa and Vasconce- 
los 3^ calculated a time-varying Hurst exponent for the 
Ibovespa returns, computed in three-year time windows 
for the period 1968-2001. A similar analysis but with 
two-year time windows is shown in Fig. 1111 One sees 
from this flgure that during the 1970's and 1980's the 
curve H{t) stays well above 1/2, the only exception to 
this trend occurring around the year 1986 when H dips 
momentarily towards 1/2 — an effect caused by the launch 
of the Cruzado economic Plan in February 1986 js^- In 
the early 1990's, after the launching of the CoUor Plan, 
we observe a dramatic decline in the curve H(t) towards 
1/2, after which it remained (within some fluctuation) 
around 1/2. This fact has led Costa and Vasconcelos 
|35j | to conclude that the opening and consequent mod- 
ernization of the Brazilian economy that begun with the 
Collor Plan resulted in a more efflcient stock market, in 
the sense that H ~ 0.5 after 1990. In Fig.^] one clearly 
sees that after the launching of a new major economic 
plan, such as the Cruzado Plan in 1986 and the Collor 
Plan in 1990, the Hurst exponent decreases. This effect 
has, of course, a simple economic interpretation: a Gov- 
ernment intervention on the market is usually designed 
to introduce "anti-persistent effects," which in turn leads 
to a momentary reduction of H. Note also that only after 
1990 does the curves H{t) goes below 1/2. This finding 
confirms the scenario described above that more devel- 
oped markets (as Brazil became after 1990) tend to have 
H^l/2. 
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D. Option pricing under the FBM assumption 



where 



We have seen above that often times asset prices have 
a Hurst exponent different from 1/2. In such cases, the 
standard Black-Scholes model does not apply, since it 
assumes that returns follow a Brownian motion (H = 
0.5). A more appropriate model for the return dynamics 
would be the fractional Brownian motion (FBM) . Indeed, 
a 'fractional Black-Scholes model' has been formulated, 
in which it is assumed that the stock price S follows a 
geometric fractional Brownian motion given by 



dS = fiSdt + aSdWH, 



(151) 



where W//(t) is the standard FBM. The fractional 
stochastic differential equation above is shorthand for the 
integral equation 

S{t) = S{0) + n f S{t')dt' + cr I S{t')dWH{t'). (152) 

To make mathematical sense of this equation is, of course, 
necessary to define stochastic integrals with respect to 
Wnit)- Here it suffices to say that a fractional Ito calcu- 
lus can indeed be rigorously defined f4^ that shares (in 
an appropriate sense) many of the properties of the usual 
Ito calculus. In the context of this fractional Ito calculus 
is possible to prove that the solution to (|151(l is given by 



Sit) ^ S'(O) exp \ ^it - -aH 



1 



2,2H 



'yWHit) 



(153) 



Compare this expression with H53|l . 

One can show that the fractional Black-Scholes 
model is complete and arbitrage- free. To price deriva- 
tives with this model, one can apply the same A-hedging 
argument used before, which now leads to the 'fractional 
Black-Scholes equation. ' The result is summarized in the 
following theorem. 

Theorem 4 In the fractional Black-Scholes model, the 
price F(S, t) of a European contingent claim with pay- 
off ^{S{T)) is given by the solution to the following 
boundary-value problem 



rS 



dF 



dS 
F{T,S) 



rF = 0, 



(154) 
(155) 



[Compare with 

For the case of a European call option the solution to 
the problem above can be found in closed form. Alterna- 
tively, one can obtain the option pricing formula directly 
from the equivalent martingale measure, without having 
to solve the above PDE. The final result for this 'frac- 
tional Black-Scholes formula' is given below 01 ■ 

Theorem 5 In the fractional Black-Scholes model, the 
price of a European call option with strike price K and 
maturity T is given by 



C{S,t) = SN{di) - Ke-''^'^-'^N{d2), 



di 



do = 



ln(f) +r(r-t) + l(7^(r^^-t^^) 

\n{^)+r{T-t)-^a^T^" -t^") 
a^/T^" - t^" 



(157) 
(158) 
(159) 



[Compare with 

The fractional Black-Scholes formula has been applied 
to price some options traded on the Brazilian market ■ 
Here, however, the option prices obtained with the for- 
mula above resulted considerably higher than those from 
the usual Black-Scholes formula. The practical relevance 
of the fractional Black-Scholes model to real markets thus 
needs to be investigated further. 



E. Multifractality in Finance 

The fact that the Hurst exponents of financial data of- 
ten display considerable variability in time indicates, as 
already mentioned, that such time series cannot be satis- 
factorily modeled in terms of a fractional Brownian mo- 
tion, which is characterized by a constant H and would 
thus capture only a sort of average behavior of the ac- 
tual price dynamics [35]. In such cases, it would be more 
appropriate to model the data as a multifractal process. 

The notion of a multifractal was first introduced in the 
context of dynamical systems to describe physical pro- 
cesses taking place on a fractal support [43 ■ A rigorous 
exposition of multifractal measures is beyond the scope 
of the present notes, and so we will content ourselves with 
a rather intuitive description of multifractality. The ba- 
sic idea here is that a monofractal process, such as the 
FBM, is characterized by a single exponent H, whereas 
for a multifractal a whole family (spectrum) of exponents 
is necessary, one for each moment of the distribution. 

We have seen above that the FBM is a Gaussian pro- 
cess whose standard deviation scales with time as 



E[W^^it)]^t". 



(160) 



If we now introduce the generalized Hurst exponents Hq 
as the corresponding scaling exponent for the 2q-th mo- 
ment, that is, 



{E[wpit)]y^'' =cqt"^, 



(161) 



where Cq is a constant, it then immediately follows from 
property H19|l of the Gaussian distribution that 



Hq^H. 



(162) 



(156) 



That is, all higher-order Hurst exponents of the FBM 
are equal to H itself and hence the FBM is said to be a 
monofractal. Our working definition of a multifractal will 
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the original literature. In turbulence |43|, there is an en- 
3- - ergy cascade from the large scales (where energy is fed 

into the system) down to the small scales (where energy 
J, o ° " is dissipated by viscosity). In currency markets the ana- 

o " log would be a kind of information cascade in time from 

^ „ ° ° long-term investors (that lake a longer view of the mar- 

,y' " ket) to short-term investors (that watch the market very 

,4° frequently) (llj . 



VIII. CONCLUSIONS 



FIG. 12: The generalized Hurst exponents Hq as a function 
of q for the Ibovespa. The dashed line indicates the linear 
behavior of a monofractal FBM with H — 0.6. 



then be a process for which the generalized Hurst expo- 
nents Hq vary with g, or alternatively, that the quantity 
qHq does not scale linearly with q. 

In general, any method used to calculate the Hurst 
exponent H (see Sec. I VII XI can be adapted to obtain the 
generalized exponents Hq. For example, the multifractal 
generalization of the DFA consists in calculating the qth- 
order fluctuation function -F'g(T), 

In complete analogy with (|148|l . the exponents Hq are 
then obtained from the scaling 

FqiT)^CqT"^. (164) 

[We remark parenthetically that the multifractal DFA 
defined above is slightly different from the formulation 
introduced in Ref. 46], but such minor distinctions do 
not matter in practice.] 

As an illustration of the multifractal DFA, we have ap- 
plied this method to the Ibovespa returns. For each q we 
computed the function Fq{T) defined in H164|l . plotted it 
in a double-logarithmic scale, and obtained the general- 
ized exponent Hq from the slope of the curve. In Fig. 1121 
it is plotted the resulting quantity qHq as a function of q. 
In this figure we clearly see that qHq deviates from the 
linear behavior expected for a monofractal, thus indicat- 
ing that the time series of Ibovespa returns does indeed 
display multifractal behavior. Evidences of multifractal 
behavior have been seen in several other stock indexes 

m. 

Multifractality has also been observed in many time 
series for currency exchange rates, and this has moti- 
vated the suggestion that there might perhaps be a for- 
mal analogy between turbulent flows and foreign cur- 
rency exchange markets; see Ref. for references to 



In these notes, I tried to present a basic introduction to 
an interdisciplinary area that has become known, at least 
among physicists working on the held, as Econophysics. 
I started out by giving some basic notions about financial 
derivatives and illustrated how to price them with a sim- 
ple binomial model. After this motivational introduction, 
I offered a concise description of Brownian motion and 
stochastic calculus, which provide the necessary mathe- 
matical tools to describe financial asset prices in continu- 
ous time. I then proceeded to discuss the Standard Model 
of Finance (SMF), namely, the Black-Scholes model for 
pricing financial derivatives. The formulation of the Ef- 
ficient Market Hypothesis, which lies at the heart of the 
SMF, in terms of martingales and its consequences for 
pricing derivatives were also discussed. Finally, I briefly 
reviewed some recent work done mostly, but not exclu- 
sively, by physicists that have produced evidences that 
the SMF may not fully describe real markets. In this 
context, some possible extensions of the Black-Scholes 
model were considered. 

I should like to conclude by mentioning that other al- 
ternatives approaches to the problem of pricing financial 
derivatives have been proposed by physicists, using meth- 
ods originally developed to treat physical problems. For 
instance, the option pricing problem was recently dis- 
cussed in the context of the so-called non-extensive sta- 
tistical mechanics ji^ . A "Hamiltonian formulation" for 
this problem was also given in which the resulting "gener- 
alized Black-Scholes" equation is formally solved in terms 
of path integrals . We refer the reader to Ref. for 
a brief review of these and other recent developments in 
Econophysics. 
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APPENDIX A: SOME BASIC DEFINITIONS 
FROM PROBABILITY THEORY 

1. Probability space 

In probability theory one usually imagines performing 

an experiment in which chance intervenes. The occur- 
rence or nonoccurrence of such an experiment is called 
an outcome u). The set of all possible elementary out- 
comes is denoted by J7. 

An event ^ is a set of outcomes, i.e., a subset of fl. 
We are interested in attributing a probability to events 
in f2. If is finite or countable, we could introduce a 
probability P(w) for each individual outcome w and then 
define the probabihty P{A) of an event A as the sum of 
the probabilities of all outcomes that make up the event 
A: 

P{A) = J2Pi^)- 
uieA 

This procedure, however, will not work when fl is un- 
countable, i.e., f7 is a continuous space such as R, since 
in this case the probability of any particular outcome is 
zero. Furthermore, a typical event will have uncountably 
many (i.e. a continuum of) outcomes. Hence the formula 
above is not applicable. That's why we need the notion 
of a probability measure to be defined shortly. 

To do this, first we need to specify the class of 'observ- 
able events', i.e., the subsets of fl to which a probability 
can be associated. If O is finite or countable, a partition 
of f2 would be the natural candidate. (Recall that a par- 
tition {Ai} of a set O is a collection of disjoint subsets, 
i.e., Ai C ^ and Ai H Aj =0 for i ^ j, whose union 
covers the whole set fi, i.e., [J^ Ai = f2.) In the case of a 
continuous space this is not possible and a different class 
of subsets is in order. To be useful, such a class must be 
closed under the various set operations, such as union, 
intersection, complementarity, etc. This is done through 
the concept of a a-algebra. 

Definition 10 A family T of subsets ofil is a a-algebra 

on CI if the following conditions are fulfilled: 

1. DeJ^andneJ^ 

2. A ^ T ^ A" ^ T, where A" = Q.\A is the 
complement of A in fl 

oo 

3. Ai,A2,...&T ^\jAi&T 

i=l 

The par (fi,^) is called a measurable space. 

The elements A C T oi the cr-algcbra are called 
measurable sets or simply events. The idea here is that 
for a given set A e ^ it is possible to ascertain whether 
any outcome u; belongs or not to A, and in this sense the 
event A is observable. 



The smallest a-algebra consists of the empty set and 
the set V, itself, i.e., J-min = {0,^^}- The largest cr- 
algebra, on the other hand, is made up of all subsets 
of f2, which is known as the power set of Q and denoted 
by 2^, hence ^max = 2^. Intermediate cr-algebras can 
be generated in the following way. Start with a given 
family U of subsets of f2 and form the intersection of all 
(7-algebras that contain U: 

= f]{:F I J- D Z^}. (Al) 

In other words, is the smallest algebra that contains 
U as is called the algebra generated by U. 

We can now attribute a 'probability of occurrence' to 
events A G via a probability measure on {Cl,J^). 

Definition 11 ^ probability measure P on the measur- 
able space (fi, JF) is a function P : —^ [0,1] such that 

1. P(0) = and P{n) = 1 

2. If Al, A2, ... £ T is a disjoint collection of elements 
ofT, i.e., AiUAj =^ifi^j, then P{U°^iAi) = 

The triple (O, J^, P) is called a probability space. 

2. Random variables 

Intuitively, a random variable AT is a function that 
attributes to each outcome uj a real number x, i.e., 
X{co) = X. We usually think of A as a random num- 
ber whose value is determined by the outcome u). A 
more formal definition is given in terms of measurable 
functions. 

Definition 12 Let {Vl,T,P) be a probability space. A 
function f : CI ^ TL is measurable with respect to the 
a-algebra T , or more compactly, T -measurable, if 

f-\u) = {Loe n\fiuj) e [/} G ^, 

for all open sets C/ € R. 

The definition above means that for any given interval 
(a, 6) C R there is a meaningful event A in fl. In an 
abuse of language we usually refer to this event as A = 
{a<f< b}. 

Definition 13 ^ random variable X on a probability 

space (f2, P) is a T -measurable function. 

A random variable X naturally generates a cr-algebra. 
This is the algebra generated by all the sets X~^{U), 
J7 C R open, and is denoted Tx- We think of Tx as 
representing the 'information' generated by the random 
variable X. The cr-algebra J^x contains the essential in- 
formation about the structure of the random variable X; 
it contains all sets of the form {uj\a < X{co) < b}. 



27 



Wc also recall the definition of the probability distri- 
bution F(x) of a random variable X: 

F{x) = P{X <x), for X e R. (A2) 

Random variables can be either discrete, if the only as- 
sume a finite or countably number of values xi,X2,---, 
or continuous. Most continuous distributions of interest 
have a density f{x), i.e., a function f{x) such that 

F{x) = f f{x)dx. (A3) 

J —CO 

This allows us to compute the probability of a given event 
A= {a < X < b} explicitly through the formula 

P{A)= f f{x)dx. (A4) 

J a 

3. Stochastic processes 

Intuitively, a stochastic process represents a dynamical 
system which evolve probabilistically in time. A formal 
definition is given below. 

Definition 14 A stochastic process is a collection of ran- 
dom variables 

defined on some probability space {Cl,J^,P) and 
parametrized by the variable t. 

We usually think of the label t as being time, so that Xt 
would represent the (random) value of the quantity X, 
say, the price of a risky asset or the position of a Brownian 
particle, at time t. For most of the cases, we consider T to 
be the halfline [0, oo). In this case we have a continuous- 
time process. Eventually, we shall also consider discrete- 
time processes, in which case the variable t assumes (non- 
negative) integer values, i.e., T = {0, 1, 2, ...}. 



It is perhaps worth emphasizing that a stochastic pro- 
cess is a function of two variables. For a fixed time i, Xt 
is a function of the random variable u, i.e., Xt — Xt{u). 
For a fixed outcome w € J7, it is a function of time, 
Xt ~ Xt{Ltj), t G T. This function of time is called a 
realization, path, or trajectory of the stochastic process 
Xt- Note, in particular, that in this context an out- 
come corresponds an entire realization or trajectory of 
the stochastic process X. 

A stochastic process is usually described in terms of the 
distributions it induces. The finite-dimensional distribu- 
tions of the stochastic process Xt are the joint probability 
distributions p{xi,ti; ...,Xn,tn) of the random variables 
Xt^, ...,Xt^, for all possible choices of times ti, ...,tn € T 
and every n > 1. The finite-dimensional distributions 
determine many (but not all) relevant properties of a 
stochastic process. A special class of stochastic processes 
are the stationary ones. 

Definition 15 A stochastic process is said to be station- 
ary if all its finite- dimensional distributions are invariant 
under a time translation, that is, 



p{xi,ti -\- T; ...,Xn,tn + t) = p(xi,ti; ...,Xn,tn), 



for any t > 0. 



Another important class of stochastic process are 
Gaussian processes, where all finite-dimensional distri- 
butions are (multivariate) Gaussians. 
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